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Abstract 

This review describes the role of magnetic symmetry in 2+1 dimen¬ 
sional gauge theories. In confining theories without matter fields in fun¬ 
damental representation the magnetic symmetry is spontaneously broken. 
Under some mild assumptions, the low-energy dynamics is determined 
universally by this spontaneous breaking phenomenon. The degrees of 
freedom in the effective theory are magnetic vortices. Their role in con¬ 
fining dynamics is similar to that played by pions and a in the chiral 
symmetry breaking dynamics. I give an explicit derivation of the effec¬ 
tive theory in (2+l)-dimensional weakly coupled conhning models and 
argue that it remains qualitatively the same in strongly coupled (2+1)- 
dimensional gluodynamics. Conhnement in this effective theory is a very 
simple classical statement about the long range interaction between topo¬ 
logical solitons, which follows (as a result of a simple direct classical cal¬ 
culation) from the structure of the effective Lagrangian. I show that if 
fundamentally charged dynamical fields are present the magnetic symme¬ 
try becomes local rather than global. The modifications to the effective 
low energy description in the case of heavy dynamical fundamental mat¬ 
ter are discussed. This effective lagrangian naturally yields a bag like 
description of baryonic excitations. I also discuss the fate of the magnetic 
symmetry in gauge theories with the Chern-Simons term. 
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1 Introduction 

This review is devoted to detailed discussion of the magnetic Zm symmetry 
in 2+1 dimensional SU{N) gauge theories. The main motivation behind 
it, is understanding confinement in gluodynamics in terms of universal 
properties determined by realization of global symmetries. 

Understanding of the low energy structure of strongly interacting the¬ 
ories is a very difficult matter. There are instances however when we 
can understand main features of the low-energy dynamics without being 
able to solve the “microscopic” theory in detail. This happens whenever 
we are lucky enough to have a spontaneously broken global continuous 
symmetry. The Goldstone theorem assures us that such a spontaneous 
symmetry breaking results in the occurrence of massless particles. The 
latter are, of course, natural low-energy degrees of freedom. Moreover, the 
original symmetry must be manifest in the effective theory that governs 
the interactions between the Goldstone bosons, albeit its implementation 
is nonlinear. This severely constrains the interactions of the Goldstone 
bosons and, in fact, lends considerable predictive power to the effective 
theory. 

The classic example of this type is the spontaneous chiral symmetry 
breaking which results in the the low-energy physics being dominated by 
pions which are described by an effective chiral Lagrangian. In real life 
the chiral symmetry is broken explicitly, the pions are massive. However, 
the effects due to the small explicit breaking are easily accommodated 
in the chiral Lagrangian. The great thing about the chiral Lagrangian, 
of course, is that it is universal. It does not care what exact dynamics is 
responsible for the symmetry breaking, what are the degrees of freedom of 
the “microscopic” theory or any other fine details. All you need to know is 
that there was a symmetry, and this symmetry was spontaneously broken. 

Can some universal considerations of a similar kind determine the 
structure of low-energy theory in pure gluodynamics? The chiral symme¬ 
try is of no relevance here. Nevertheless, I will argue that the symmetry 
path is a very fruitful one. The main thesis is this: there is a global dis¬ 
crete symmetry in gluodynamics, which is spontaneously broken in the 
vacuum. In what follows this symmetry will be called the magnetic Zn 
symmetry. Although it is discreet and, therefore, does not have all the 
bliss of the Goldstone theorem, under some natural assumptions it does 
indeed determine the low-energy dynamics. In physical (3+l)-dimensional 
theory this symmetry is of a somewhat unusual type - its charge is not 
a volume integral, but rather a surface integral. The implications of such 
a symmetry for the low-energy Lagrangian have to be studied in more 
detail; this has not been done so far. However, one can go much further 
in 2+1 dimensions, where this symmetry is of a familiar garden variety. 
This review will, therefore, be devoted to (2+l)-dimensional non-Abelian 
theories. The structure of the paper is the following. In Sec. 2 I will 
discuss in detail the nature of the magnetic Zn symmetry and the ex¬ 
plicit construction of both, the generator of the group and the local order 
parameter. In Sect. 3 I will show how this symmetry is realized in dif- 
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ferent phases of Abelian gauge theories. I will show that in the Coulomb 
phase this symmetry is spontaneously broken and in the Higgs one it is 
respected by the vacuum state. In the Abelian theories the magnetic sym¬ 
metry is presented by a continuous U(l) group - its spontaneous breaking 
leads to the appearance of a massless particle, the photon. I will derive 
the low-energy effective Lagrangian that describes this symmetry breaking 
pattern and show that it exhibits logarithmic confinement in the Coulomb 
phase. In Sec. 4 I will discuss non Abelian confining theories. Here the 
magnetic symmetry is discrete, but it is also spontaneously broken as 
in the Coulomb phase of QED. In weakly interacting theories, like the 
Georgi-Glashow model, it is possible to derive an effective Lagrangian 
which follows from this symmetry breaking pattern. I will demonstrate 
that, due to the fact that the magnetic symmetry is discrete, the effective 
theory exhibits linear confinement, and that this confinement mechanism 
in the effective theory appears very simply at the classical level. I will give 
arguments that the basic structure of the effective Lagrangian as well as 
physics of confinement stay the same in the limit of pure gluodynamics. 
I will discuss similarities and differences between the confining properties 
of the weakly interacting non-Abelian theories and gluodynamics. In Sec. 
5 I will discuss the modifications to the magnetic symmetry concept when 
fundamentally charged matter fields are present. It turns out that in the¬ 
ories of this kind the magnetic Zn symmetry is gauged. Consequently the 
low energy effective theories is a local Zn invariant model. It is however 
still quite useful. In particular this approach yields a distinctive low en¬ 
ergy description of baryonic states as ’’bags” of deconfined phase which 
contain the fundamental quarks inside them. Sec. 6 is devoted to the 
discussion of SU{N) gauge theories with the Chern-Simons term. I show 
that the magnetic symmetry is unbroken in these theories in the pertur¬ 
bative domain. It is however possible that it is spontaneously broken 
at low values of the Chern-Simons coefficient, and the theories are then 
confining. I end with a short discussion in Sec.7. 


2 The magnetic symmetry in 2+1 di¬ 
mensions 

A while ago ’t Hooft gave an argument establishing that a non-Abelian 
SU(A) gauge theory without charged fields in the fundamental represen¬ 
tation possesses a global Zn symmetry. [Q] 

Consider first a theory with several adjoint Higgs fields, so that by 
varying parameters in the Higgs sector the SU(Af) gauge symmetry can 
be broken completely. In this phase the perturbative spectrum will contain 
the usual massive “gluons” and the Higgs particles. However, in addition 
to those, there will be heavy stable magnetic vortices. Those are analogs 
of Abrikosov-Nielsen-Olesen vortices in the superconductors; they must 
be stable by virtue of the following topological argument. All the fields 
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in the vortex configuration away from the vortex core are pure gauge, 

H°‘{x) = U{x)h", A^ = iUd>^UK (1) 

Here the index a labels the scalar fields in the theory, are the constant 
vacuum expectation values of these fields, and U(x) is a unitary matrix. 
As one goes around the location of the vortex in space, the matrix U winds 
nontrivially in the gauge group. This is possible, since the gauge group in 
the theory without fundamental fields is S\]{N)/Zn- It has a nonvanishing 
first homotopy group 7ri(SU(A^)/Z]v) = Zn- Practically, this means that 
when one makes a full circle around the vortex location, U does not return 
to the same SU(A'^) group element C/q, but, rather, ends up at exp{i^}[/o- 
Adjoint fields do not feel this type of discontinuity in U and, therefore, 
the energy of such a configuration is finite. Since such a configuration 
can not be smoothly deformed into a trivial one, a single vortex is stable. 
Processes involving annihilation of N such vortices in the vacuum are 
allowed since the A^-vortex configurations are topologically trivial. One 
can, of course, find explicit vortex solutions once the Higgs potential is 
specified. As any other semiclassical solution at weak coupling, its energy 
is inversely proportional to the gauge coupling constant and is, therefore, 
very large. One is faced with the situation where the spectrum of the 
theory contains a stable particle even though its mass is much larger than 
the masses of many other particles (the gauge and Higgs bosons) and 
the phase space for its decay into these particles is enormous. The only 
possible reason for the existence of such a heavy stable particle is that it 
carries a conserved quantum number. The theory must, therefore, possess 
a global symmetry which is unbroken in the completely Higgsed phase. 
The symmetry group must be Zjv since the number of vortices is only 
conserved modulo N. ’t Hooft dubbed this symmetry “topological.” I 
prefer to call it “magnetic” for reasons that will become apparent in a 
short while. 

Now imagine smoothly changing the parameters in the Higgs sector 
so that the expectation values of the Higgs fields become smaller, and, 
finally, the theory undergoes a phase transition into the confining phase. 
One can further change the parameters so that the adjoint scalars become 
heavy and eventually decouple completely from the glue. This limiting 
process does not change topology of the gauge group and, therefore, does 
not change the symmetry content of the theory. One concludes that the 
pure Yang-Mills theory also possesses a Zn symmetry. Of course since the 
confining phase is separated from the completely Higgsed phase by a phase 
transition one may expect that the Zn symmetry in the confining phase 
is realized differently in the vacuum than in the completely “Higgsed” 
phase. In fact, the original paper by ’t Hooft, as well as subsequent 
work,Q convincingly argued that in the confining phase the Zn symmetry 
is spontaneously broken, the breaking being related to the confinement 
phenomenon. We will have more to say about this later. 

The physical considerations given above can be put on a firmer formal 
basis. In particular, one can explicitly construct [B the Zn generator as 
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well as the order parameter associated with it, the operator that creates 
the magnetic vortex. We start this discussion by considering Abelian 
theories, where things are simpler and are under full control. 


2.1 Abelian theories 

In the U(l) case the homotopy group is Z and, therefore, the magnetic 
symmetry is U(l) rather than Zn- Identifying the relevant charge is, in 
fact, an absolutely straightforward task. It is nothing but the magnetic 
flux through the equal time plane, with the associated conserved current 
being the dual of the electromagnetic field strength 

-I- = / cPxB{x ), d'^F^ = 0. (2) 


The current conservation is insured by the Bianchi identity. It may 
come as a surprise that we are seriously considering a current whose con¬ 
servation equation is an “identity.” However “identity” is a relative thing. 
The conservation equation is trivial only because we have written the com¬ 
ponents of the field strength tensor in terms of the vector potential. The 
introduction of the vector potential is none other than the potentiation of 
the conserved current Ffj ,, that is explicit solution of the conservation equa¬ 
tion. In exactly the same way one can potentiate any conserved current, 
and such a potentiation will turn the pertinent conservation equation into 
identity. Thus has exactly the same status as any other local conserved 
gauge invariant current, and should be treated as such. 

Once we have the current and the charge, we also know the elements 
of the symmetry group. A group element of the U(l) magnetic symmetry 
group is 

Wa(oo) = exp{iQ;$} (3) 

for the arbitrary value of a. The notation W is not accidental here, 
since the group element is indeed a large spatial Wilson loop defined on 
a contour that encloses the whole system. 

The question that might bother us is whether this group acts nontriv- 
ially on any of the local physical observables in our theory. The obvious 
gauge invariant observables like B and E commute with IT. However, 
there is another set of local gauge invariant observables in the theory 
which do indeed transform nontrivially under the action of IT. Consider, 
following ’t Hooft, the operator of the “singular gauge transformation,” 


V{x) = exp ■ 


j2 

d y 


jx - y)j 

{x - yY 


Ei{y) + Q{x - y),h{y) 


(4) 


where Q{x — y) is the polar angle function and Jo is the electric charge 
density of whatever matter fields are present in the theory. The cut discon¬ 
tinuity in the function 0 looks bothersome; but it is, in fact, not physical 
and completely harmless. The gauge function jumps across the discon¬ 
tinuity by 27r/gr, but since the only dynamical fields in the theory have 
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charges that are integer multiples of g, the discontinuity is not observable. 
The cut can be chosen parallel to the horizontal axis. Using the Gauss’ 
law constraint this can be cast in a different form, which we find more 
convenient for our discussion, 


V(x) = exp — f dy^ e-ijEi{y ), (5) 

9 Jc 

where the integration goes along the cut of the function 0 which starts 
at the point x and goes to spatial infinity. In this form it is clear that the 
operator does not depend on where precisely one chooses the cut to lie. 
To see this, note that changing the position of the cut C to C adds 

— / d^xdiE' 

9 Js 

to the phase (here S is the area bounded by C — G'). In the theory 
under consideration only charged particles where charges multiples of g 
are present. Therefore, the charge within any closed area is a multiple 
integer of the gauge coupling d^xdiE^ = gn and the extra phase factor 
is always unity. The only point in space where the action of V{x) on any 
physical state is nontrivial is the point x. Therefore, the field V {x) acts 
like any other local field. With a little more work one can prove not 
only that V{x) is a local field, in the sense that it commutes with any 
other local gauge invariant operators 0(y), x ^ y, but also that it is a 
bona fide Lorentz scalar. 

The physical meaning of the operator V is very simple. Calculating 
its commutator with the magnetic field B we find 

Vix)Biy)V\x) = B{y) + ^5^{x-y). (6) 

Thus, V creates a pointlike magnetic vortex of flux ‘1-k jg. This commuta¬ 
tor also tells us that 


WlV{x)W^ = e"^V{x ), (7) 

and, therefore, V is indeed a local eigenoperator of the magnetic U(l) 
symmetry group. 

Equations (g) and (|^ formalize the physical arguments of’t Hooft in 
the Abelian case. We have explicitly constructed the generator and the 
local order parameter of the magnetic symmetry. 

2.2 Non-Abelian theories at weak coupling 

Let us now move on and consider an analogous construction for the non- 
Abelian theories. Ultimately we are interested in pure Yang-Mills theory. 
However, it is quite illuminating to start with the theory with an adjoint 
Higgs field and take the decoupling limit later. For simplicity we discuss 
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the SU(2) gauge theory. Consider the Georgi-Glashow model - SU(2) 
gauge theory with the adjoint Higgs field, 

-C = , (8) 

where 

VfH'’ = . (9) 

At large and positive the model is weakly coupled. The SU(2) gauge 
symmetry is broken down to U(l) and the Higgs mechanism takes place. 
Two gauge bosons, , acquire a mass, while the third one, the “photon,” 
remains massless to all orders in perturbation theory. The theory in this 
region of the parameter space very much resembles electrodynamics of 
the vector charged field. The Abelian construction can, therefore, be 
repeated. The SU(2) gauge invariant analog of the conserved dual field 
strength is 

F'' = , (10) 

^ 9 

where is the unit vector in the direction of the Higgs field. 

Classically this current satisfies the conservation equation 

= 0 . ( 11 ) 

The easiest way to see this is to choose a unitary gauge of the form 
H°'{x) = H (a:)5“®. In this gauge F is equal to the Abelian part of the dual 
field strength in the third direction in the color space. Its conservation 
then follows from the Bianchi identity. 

Thus, classically the theory has a conserved U(l) magnetic charge 
^ ~ f d^xFo, just like QED. However the unitary gauge can not be im¬ 
posed at the points where H vanishes, which necessarily happens in the 
core of the ’t Hooft-Polyakov monopole. It is well-known, of course,^ 
that the monopoles are the most important nonperturbative configura¬ 
tions in this model. Their presence leads to a nonvanishing small photon 
mass, as well as confinement of the charged gauge bosons, with a tiny 
nonperturbative string tension. As far as the monopole effects on the 
magnetic flux are concerned, their presence leads to a quantum anomaly 
in the conservation equation ([ll|). As a result, only the discrete Z 2 sub¬ 
group of the transformation group generated by remains unbroken in 
quantum theory. The detailed discussion of this anomaly, the residual 
Z 2 symmetry and their relation to monopoles is given in Ref. 15. The 
nonanomalous Z 2 magnetic symmetry transformation is generated by the 
operator 

C/ = exp{i|$}. (12) 

The order parameter for the magnetic Z 2 symmetry is constructed analo¬ 
gously to QED as a singular gauge transformation generated by the gauge 
invariant electric charge operator, 

J d^xJo{x). 




(13) 


Explicitly 


V(x) 


i 

exp — 
9 


j2 

d y 


€ij 


jx - y)j 

{x-yY 


n"{y)Ei{y) +e(a; 


y)Jo{y) 


exp— [ dy'- eijU^-Ejiy). 
9 J c 


( 14 ) 


One can think of it as of a singular SU(2) gauge transformation with the 
field dependent gauge function 


A“(i/) = ie(a;-j/)n“(i7). 


(15) 


This field dependence of the gauge function ensures the gaMe invariance 
of the operator V. Just like in QED, it can be shown M, M that the 
operator E is a local scalar field. Again like in QED, the vortex operator 
E is a local eigenoperator of the Abelian magnetic field B{x) = Fq. 

97r 

[V{x),B{y)] = - — V{x)5\x-y). (16) 

That is to say, when acting on a state it creates a pointlike magnetic 
vortex which carries a quantized unit of the magnetic flux. 

The magnetic Z 2 acts on the vortex field V as a phase rotation by tt 


e'^^V{x)e- 


= -V(x). 


(17) 


This is an explicit realization of the magnetic Z 2 symmetry in the Georgi- 
Glashow model. 


2.3 Gluodynamics 

From the Georgi-Glashow model we can easily get to pure Yang-Mills 
theory. This is achieved by smoothly varying the coefficient in the La- 
grangian, so that it becomes negative and eventually arbitrarily large. In 
this limit the Higgs field has a large mass and, therefore, decouples leav¬ 
ing pure gluodynamics behind. It is well known that in this model the 
weakly coupled Higgs regime and the strongly coupled confining regime 
are not separated by a phase transition. The limit of pure Yang-Mills in 
this model is, therefore, smooth. 

In the limit of pure Yang-Mills the expressions (p^, (|l^, and 
have to be taken with care.|12|] When the mass of the Higgs field is very 
large, the configurations that dominate the path integral of the theory are 
those with a very small value of the modulus of the Higgs field \H\ oc 1/M. 
The modulus of the Higgs field, in turn, controls fluctuations of the unit 
vector n“, since the kinetic term for n in the Lagrangian is |Hp(D^n)^. 
Thus, as the mass of the Higgs field increases, the fluctuations of n grow 
both, in amplitude and frequency, and the magnetic field operator B as 
defined in Eq. (^) fluctuates wildly. Of course, this situation is not 
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unusual. It happens whenever one wants to consider (in the effective low- 
energy theory) an operator which explicitly depends on fast, high-energy 
variables. The standard way of dealing with it is to integrate over the fast 
variables. There could be two possible outcomes of this “integrating out” 
procedure. Either the operator in question becomes trivial (if it depends 
strongly on the fast variables), or its reduced version is well-defined and 
regular in the low-energy Hilbert space. The “magnetic field” operator B 
in Eq. (|^ is obviously of the first type. Since in the limit of pure Yang- 
Mills all orientations of n“ are equally probable, integrating over the Higgs 
field at fixed will lead to vanishing B. However, what interests us is 
not so much the magnetic field but, rather, the generator of the magnetic 
Z 2 transformation. It is actually instructive to consider an operator that 
performs the Z 2 transformation not everywhere in space, but only inside 
a contour C 

Uc = exp{i^ / d^xB{x)}, (18) 

with the area S being bounded by C. In the limit of gluodynamics we are 
led to consider the operator 


Uc 


lim 

H^O 


Dn°‘ exp I — \H\^{Dna)^'^ 


X exp|i^ J — ie*-’e“*’“na(T’in)*’(r>jn)“) I .(19) 

The weight for the n integration is the kinetic term for the isovector n^. 
As was noted before, the action does not depend on n“ in the Yang-Mills 
limit since —> 0. However, the first term in Eq. regulates the path 
integral and we keep it for this reason. This operator may look somewhat 
unfamiliar, at first sight. Hoover, in a remarkable paper (see Ref. 18) 
Diakonov and Petrov showed □ that Eq. © is equal to the trace of the 
fundamental Wilson loop along the contour C. 

Uc = Wc = TrP exp jig J . (20) 


Taking the contour C to run at infinity, we see that in gluodynamics 
the generator of the magnetic Z 2 symmetry is the fondamental spatial 
Wilson loop along the boundary of the spatial plane.u 

^ We note that Diakonov and Petrov had to introduce a regulator to define thp na.th integral 
over n. The regulator they required was of precisely the same form as in Eq. (|l9|) . 

^There is a slight subtlety here that may be worth mentioning. The generator of a unitary 
transformation should be a unitary operator. On the other hand, the trace of the fundamental 
Wilson loop is not unitary. Therefore, one should strictly speaking consider, instead, a unita- 
rized Wilson loop W = , ^ However, the factor between the two operators VWW^ is an 

operator that is only sensitive to the field behavior at infinity. It commutes with all physical 
local operators 0{x) unless x —s- 00 . In this aspect, it is very different from the Wilson loop 
itself, which has a nontrivial commutator with vortex operators V{x) at all values of x. Since 
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Of course, one does not have to go through the exercise with the 
Georgi-Glashow model in order to show that the fundamental Wilson loop 
generates a symmetry. Instead one can directly consider the commutator 


[W, H] = lim 

C —►oc 


0 . ( 21 ) 


dxiTrVEi{x)exp{i (b dj/iA*(i/)} —>< 

ic JC(x,x) 

Here the integral in the exponential on the right-hand side starts and 
ends at the point of insertion of the electric held. For a hnite contour C 
the commutator does not vanish only along the contour itself - it does 
not contain any bulk terms. Making the contour C go to inhnity and 
assuming, as usual, that in the theory with a hnite mass gap at inhnity 
no physical modes are excited we condude that the commutator of H with 
inhnitely large Wilson loop vanishes.0 

Next we consider the vortex operator Eq. ( 0 )- Again, in order to 
hnd the limit of pure Yang-Mills we have to integrate this expression over 
the orientations of the unit vector n“. This integration is equivalent to 
averaging over the gauge group. Following Ref. 18 one can write Ua in 
terms of the SU(2) gauge transformation matrix Q, 


n = -TrHrHVs . (22) 

The vortex operator in the limit of pure gluodynamics then becomes 


V{x) = 


DO. exp ^ 


dyitijTxQ.EjQ} Tz . 


(23) 


This form makes it explicit that V (a:) is defined as the gauge singlet part 
of the following, apparently non-gauge-invariant operator: 


V (x) = exp 


27ri 




(24) 


the correlators of all gauge invariant local fields in the pure Yang-Mills theory are massive and, 
therefore, short range, the operator y^tT(C)iyi(C) at C —»■ oo must be a constant operator 
on all finite-energy states. Therefore, the difference between W and IT is a trivial constant 
factor and we will not bother with it in what follows. 

^Note that the nonvanishing of the commutator at finite C is of exactly the same nature 
as for any other “conserved charge,” which is defined as an integral of local charge density 
Q = lim(7^oo J"| d'^xp{x). The commutator of such a charge with the Hamiltonian also 

contains surface terms, since the charge density p itself never commutes with the Hamiltonian. 
Rather, the commutator is a total derivative. For a conserved charge, due to the continuity 
equation, this surface term is equal to the circulation of the spatial component of the current 


[Q,H]=^c- 


dx'-ji 


Again, the vanishing of this term is the consequence of the vanishing of the physical fields at 
infinity in the theory with a mass gap. When the charge is not conserved, the commutator 
contains a bulk term, in addition to the surface term. It is the absence of such bulk terms that 
is a unique property of the conserved charge. The same conclusion is reached if, rather than 
considering the generator of the algebra, one considers the commutator of the group ele ment 
for either continuous or discrete symmetry groups. Therefore, the commutator in Eq. ( [21] ) 
indeed tells us that W is a conserved operator. 
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The integration over obvionsly projects out the gange singlet part of V. 

In the present case, however, this projection is redundant. This is because 
even though V itself is not gauge invariant, when acting on a physical state 
it transforms it into another physical state.u By physical states we mean 
the states which satisfy Gauss’ constraint in pure Yang-Mills theory. This 
property of V was noted by’t Hooft.0 To show this let us consider V(x) 
as defined in Eq. ( |2^ ) and its gauge transform Vh = 17^1/II where 17 is an 
arbitrary nonsingular gauge transformation operator. The wave functional 
of any physical state depends only on gauge invariant characteristics of 
the vector potential, i.e. only on the values of the Wilson loops over all 
possible contours, 

= ']?[{IE(C)}]. (25) 

Acting on this state by the operators V and Vh, respectively, we obtain 

E|T) = Tv[Ai] = T[{ytT(C7)E+}], 

Yn|v&) = (26) 

However, it is easy to see that the action of V(x) and Vh(a;) on the Wilson 
loop is identical - they both multiply it by a phase belonging to the center 
of the group if x is inside C, and do nothing otherwise. Therefore, 

Y|^) = Vh|^) (27) 

for any physical state <&. Thus, we have 

HHIT) = 17^17^1^) = HIT), (28) 

where the first equality follows from the fact that any physical state is 
invariant under the action of any gauge transformation 17, while the second 
equality follows from Eq. (p^. This equation is nothing but the assertion 
that the state H|T) is physical, i.e. invariant under any nonsingular gauge 
transformation. 

Thus, we have proved that when acting on a physical state the vortex 
operator creates another physical state. Given such an operator, the gauge 
invariant projection only affects its matrix elements between unphysical 
states. Since we are interested only in calculating correlators of V between 
physical states, the gauge projection is redundant and we can freely use 
V rather than V to represent the vortex operator. 

The formulae in this section can be straightforwardly generalized to the 
SU(A^) gauge theories. Once again one can start with the Georgi-Glashow- 
like model, where SU(A’) is Higgsed to U(l)^'^“^\ The construction of 
the vortex operator and the generator of Zn in this case is very similar; 
necessary details are given in Ref. 15. Taking the mass of the Higgs held 
to inhnity again projects the generator onto the trace of the fundamental 
Wilson loop. The vortex operator can be taken as 

H(a;) = exp{^ J dy'^eijTr{YEi{y))} (29) 

“^This is not a trivial statement, since a generic non-gauge-invariant operator has nonvan¬ 
ishing matrix elements between the physical and unphysical sectors. 
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where the hypercharge generator Y is defined as 

y = diag(l,l,...,-(iV-l)) , (30) 

and the electric field is taken in the matrix notation Ei = A“i7“. Here A“ 
stand for the SU(Ai) generator matrices in the fundamental representa¬ 
tion. 

To summarize this section, we have established two important facts. 
First, SU(Ai) gauge theories in 2-1-1 dimensions have a global Zn mag¬ 
netic symmetry. The generator of this magnetic symmetry group is the 
fundamental Wilson loop around the spatial boundary of the system. Sec¬ 
ond, this symmetry has a local order parameter. This order parameter 
is a local gauge invariant scalar field which creates a magnetic vortex of 
fundamental flux. 

The next question we should ask is whether this global symmetry is 
at all relevant for low-energy dynamics. In the next section we will show 
that this is indeed the case. We will calculate the expectation value of V 
in confining and nonconfining situations and will show that confinement 
is rigidly related to the spontaneous breaking of the magnetic symmetry. 


3 The vacuum realization of the mag¬ 
netic symmetry, the effective Lagrangian 
and confinement. Abelian theory 

Again, we start our discussion with the Abelian theory. Consider a U(l) 
gauge theory with scalar matter field 

£ = - A(^»" . (31) 

Depending on the values of the coupling constant this theory can be ei¬ 
ther in the Coulomb phase with the massless photon and logarithmically 
confined charges, or in the Higgs phase, which is massive - all electric 
charges are screened. 

3.1 Realizations of the magnetic symmetry 

Let us start by calculating (V) in the Coulomb phase. This can be done 
using the standard weak coupling perturbation theory. The expectation 
value of V is given by the following expression: 

(!/(*)) = V-i hm / 

T—*-oo / 

(32) 

Here N ^ is the normalization factor - the usual vacuum-to-vacuum am¬ 
plitude, |0) is the perturbative Fock vacuum and the integral over Aq is 
the standard representation of the projection operator which projects the 
Fock vacuum |0) onto the gauge invariant subspace which satisfies Gauss’ 
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law. As usual, discretizing time, introducing resolution of identity at ev¬ 
ery time slice and integrating over Ei in the phase space path integral, 
this expression can be rewritten as a path integral in the field space. The 
result is easy to understand - it is almost the same as for the vacuum-to- 
vacuum amplitude, except that at time t = 0 the spatial derivative of the 
scalar potential Aq is shifted by a c-number field due to the presence of 
the vortex operator. 


{V{x)) = N-^ 


VA^ exp 


1 

4 


d^y[Ft,{y) - U{y - x)f + Lmggs 


(33) 

The c-number field is the magnetic field of an infinitely thin magnetic 
vortex which terminates at the point x. One can view it as the Dirac 
string of a (three-dimensional Euclidean) magnetic monopole, 


" ~ ~ ZTT 

/o = /2 = 0 My) = — 0iyi)S{y2)S{y3). (34) 

g 

Thus, at weak coupling we have to find the solution of the classical 
equations of motion following from the action with the external source 
Eq. (|3^). The nature of this solution is clear: it is just the Dirac monopole. 
The action of this solution is infrared-finite, since the contribution of 
the Dirac string (which normally would be linearly infrared-divergent) is 
cancelled by the external source, 


{V) = exp{-Sci} , (35) 

with 

5ci = 4 • (36) 

Here A is an ultraviolet cutoff which has to be introduced since the ac¬ 
tion of the pointlike monopole diverges in the ultraviolet. This ultravio¬ 
let divergence is benign since it can be eliminated by the multiplicative 
renormalization of the vortex operator.The important point is that 
since there is no divergence in the infrared, the expectation value of V is 
nonvanishing. Thus, we conclude that in the Coulomb phase of QED the 
magnetic symmetry is spontaneously broken. 

The spontaneous breaking of a continuous symmetry must be accom¬ 
panied by a massless Goldstone particle. Indeed, in QED such a particle 
exists - this is the massless photon. The matrix element of the magnetic 
current between the vacuum and the one-photon state is 

{Q\F^\ki) ^ (37) 

where Z(Q) is the on-shell photon wave function renormalization. 

This is the standard form of the matrix element of a spontaneously 
broken current, with Z{0) playing the role of /^. 

Let us now perform the same calculation in the Higgs phase. The 
path integral representation in Eq. ( |3^ is still valid. However, the classi¬ 
cal solution that dominates this path integral is now very different. Since 
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in the Higgs phase the photon has a nonzero mass /i the classical ac¬ 
tion of the three-dimensional monopole in the superconducting medium 
is linearly divergent in the infrared. Q Essentially, the magnetic flux that 
emanates from the monopole can not spread out in space (time) but is 
rather concentrated inside a flux tube of the thickness 1/^ which starts at 
the location of the monopole and goes to a spatial boundary at infinity. 
The action of such a held conhguration is proportional to the linear size 
of the system and diverges in the thermodynamic limit. As a result, the 
expectation value of V in the Higgs phase vanishes, 

(V) = e-^ 0 ■ (38) 

Similarly, the vortex held correlator is given in terms of the classical energy 
of a monopole-antimonopole pair in the superconductor. The Euclidean 
action of this conhguration is proportional to the distance between the 
monopole and the antimonopole and, therefore, the correlator of the V 
decays exponentially, 

(H*(*)l/(y))~e-"^'^l“^-"l, (39) 

with Mv being the mass of the ANO vortex. 

This simple calculation can be improved perturbatively. In the next- 
to-leading order one has to calculate the determinant of the Schrodinger 
operator of a particle in the held of a monopole, and this corrects the value 
of the mass Mv- Higher orders in perturbation theory can be calculated 
too; we will not pursue this calculation here. 

The main lesson is that the expectation value of the order parameter 
vanishes in the Higgs phase; thus, the magnetic symmetry is unbroken. 

3.2 The low-energy effective Lagrangian and the 
logarithmic confinement of electric charges 

Thus we see that the Coulomb-Higgs phase transition can be described as 
that due to the restoration of the magnetic U(l) symmetry, the pertinent 
local order parameter being the vortex operator V. It should then be 
true that the low-energy dynamics in the vicinity of the phase transition 
is described by an effective low-energy Lagrangian. For the U(l) symme¬ 
try breaking such an effective Lagrangian can be immediately written as 
follows: 

C = d^V*d>^V-\{V*V - gf . (40) 

Although this Lagrangian may seem a bit unfamiliar in the context 
of QED, a little thought convinces one that it indeed describes all the 
relevant light degrees of freedom of the theory. In the Coulomb phase, 
where {V) = g ^ 0, the physical particles are interpolated by the phase 
and the radial part of V, 

V (x) = pe^ . (41) 

The phase x is of course the massless Goldstone boson held, i.e. the 
photon. The fluctuation of the radial component p—p is the lightest scalar 
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particle, which in this case is the lightest meson, or scalar positronium. In 
the Higgs phase, the field V interpolates physical excitations - the ANO 
vortices. Of course, far from the phase transition the vortices are heavy 
and there are other, lighter excitations in the spectrum. The validity of 
this effective Lagrangian on the Higgs phase side is therefore limited to a 
narrow critical region where vortices are indeed the lightest particles. 

The charged particles have not appeared in our discussion so far. In¬ 
deed in 2-1-1 dimensions the electrical charges are confined, and, therefore, 
we do not expect the charged fields to appear as basic degrees of freedom 
in the effective low-energy Lagrangian. However, our original purpose was 
precisely to understand the mechanism of confinement through studying 
the effective Lagrangian. Thus, if we are unable to identify the charged 
objects in this framework our program is doomed to failure. Fortunately 
it is not difficult to understand how the charged states are represented 
in the Lagrangian Eq. The easiest way to do this is to identify the 

electric charge through the Maxwell equation, 

= -^^p,v\dvF\ . (42) 

The dual field strength is obviously proportional to the conserved 
U(l) current 

= ^i{V*dV - h.c.). (43) 

The proportionality constant in this relation is dictated by the fact that 
in the Higgs phase the magnetic vortices that carry one unit of the U(l) 
charge, carry the magnetic flux of ‘I'Kjg. The above two relations give 

-Jm = jwa.(H*aAH). (44) 

TT 

To calculate the electric charge we integrate the zeroth component of the 
current over the two-dimensional plane 

— dxidiX- (45) 

Jc^oo 

Therefore, the electric charge is proportional to the winding number of 
the phase of the field V. 

So, the charged states do appear in the low-energy description in a very 
natural way. A charged state is a soliton of V with a nonzero winding 
number. This identification immediately tells us that the charged par¬ 
ticles are logarithmically confined. Consider, for example, the minimal 
energy configuration in the sector with the unit winding number. This is 
a rotationally invariant hedgehog. Fig. |^, which, far from the soliton core, 
has the form 

V{x)=ne'^^^\ (46) 

Here 9{x) is the angle between the vector x and one of the axes. 
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Figure 1: The hedgehog configuration of the field V in the state of unit charge 
in the Abelian theory. 


The self-energy of this configuration is logarithmically divergent in the 
infrared due to the contribution of the kinetic term 

E ^ . (47) 


This is nothing but the electromagnetic self-energy of an electrically charged 
state associated with the logarithmic Coulomb potential in two spatial di¬ 
mensions. Therefore, the logarithmic confinement of QED is indeed very 
easily and transparently seen at the level of the low-energy effective La- 
grangian. In itself this is, perhaps, not such a big deal, since confinement 
in this model is a kinematic phenomenon: it is a direct consequence of the 
logarithmic behavior of the Coulomb potential. We will see later, how¬ 
ever, that this low-energy picture also naturally generalizes to non-Abelian 
theories and easily accommodates linear confinement. 

Before moving on to the non-Abelian theories, let me make two com¬ 
ments. First, in perturbative regime the couplings of the effective La- 
grangian can be determined in terms of the couplings of the fundamental 
QED Lagrangian. To determine two constants in Eq. one needs two 
matching conditions. One of them can be naturally taken as the coefficient 
of the infrared logarithm in the self-energy of a charged state. Matching 
Eq. (^) with the Coulomb self-energy gives 




2 



(48) 


The other coefficient is determined by requiring that the mass of the radial 
excitation p matches the mass of the scalar positronium, which to leading 
order in g is just 2M. This condition gives 


A = 


2n^M 


(49) 


The second comment is about the relation between the vortex operator 
as defined in Eq. (^^ and the field V that enters the effective Lagrangian 
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The vortex operator, as defined in Eq. (|l4|) has a fixed length, whereas 
the field V which enters the Lagrangian (M® is a conventional complex 
field. How should one understand that? First of all, at weak gauge cou¬ 
pling the quartic coupling in the dual Lagrangian is large, A ^ oo. This 
condition freezes the radius of V dynamically. In fact, even at a finite 
value of A, if one is interested in the low-energy physics, the radial com¬ 
ponent is irrelevant as long as it is much heavier than the phase. Indeed, 
at weak gauge coupling the phase of V (which interpolates the photon) is 
much lighter than all other excitations in the theory. Effectively at low 
energies Eq. (^) reduces to a nonlinear a model and one can identify the 
field V entering Eq. (^) directly with the vortex operator of Eq. ©■ 
However, it is well known that quantum mechanically the radial degree 
of freedom of the a model field is always resurrected. The spectrum of 
such theory always contains a scalar particle which can be combined with 
the phase into a variable-length complex field. The question is purely 
quantitative - how heavy is this scalar field relative to the phase? 

Another way of expressing this is as follows. The fixed length field 
V is defined at the scale of the ultraviolet cutoff in the original theory. 
To arrive at a low-energy effective Lagrangian one has to integrate over 
all quantum fluctuations down to some much lower energy scale. In the 
process of this integrating out, the field is “renormalized” and acquires 
a dynamical radial part. Then the mass of this radial part is equal to 
the mass of the lowest particle with the same quantum numbers in the 
original theory. That is why the parameters in Eq. must be such that 
the mass of the radial part of V would be equal to the mass of the lightest 
scalar positronium. 


4 Non-Abelian theories 

In this section I want to extend the construction of the low-energy effective 
Lagrangian discussed above to the non-Abelian theories and see how the 
realization of the magnetic symmetry in this Lagrangian is related to 
confinement. Before discussing specific models in detail, let me present 
a general argument which establishes that the spontaneous breaking of 
magnetic .Zat in non-Abelian theories implies the area law behavior of the 
fundamental Wilson loop. 

4.1 Broken Zjv means confinement 

As we discussed in the Sec. 2, the generator of the magnetic Zn in pure 
gluodynamics is the fundamental Wilson loop around the spatial bound¬ 
ary of the system. By the same token, the Wilson loop around a closed 
spatial contour C generates the Zn transformation at the points inside the 
contour C. Let us imagine that the Zn symmetry is spontaneously bro¬ 
ken in the vacuum and consider the expectation value of W (C) in such a 
state |0) . The expectation value (0|IT(C)|0) is nothing but the overlap of 
the vacuum state (0| and the state |S) which is obtained by acting on the 


18 


vacuum by W{C), |S) = W{C)\Q). If the symmetry is broken, the wave 
function |0) depends explicitly on the degrees of freedom which are non¬ 
invariant under the symmetry transformation, and is peaked around some 
specific orientation of these variables in the group space. For simplicity 
let us think about all these non-invariant variables as being represented 
by the vortex field V. In the vacuum state the field V has a nonvanishing 
vacuum expectation value (VEV) and points in some fixed direction in 
the internal space. On the other hand, in the state l^) its direction in the 
internal space is different - rotated by 2-11 jN - at the points inside the 
area S bounded by C, since at these points the field V has been rotated by 
the action of W{C). In the local theory with finite correlation length the 
overlap between the two states approximately factorizes in the product 
of the overlaps taken over the regions of space of linear dimension of the 
order of the correlation length 1, 

(015) =n,(0.|S.), (50) 

where the label x is the coordinate of the point in the center of a given 
small region of space. For x outside the area S the two states |0a;) and |5a;) 
are identical and, therefore, the overlap is unity. However, for x inside 
S the states are different and the overlap is a number e~^ smaller than 
unity. The number of such regions inside the area is obviously of the order 
of S/l^ and, thus, 

(IT(C')) =exp{- 7 ^}. (51) 

Therefore, in the broken phase the spatial Wilson loop has the area law 
behavior. 

Now consider the unbroken phase. Again the average of W( C) has the 
form of the overlap of two states which factorizes as in Eq. (p0|). Now, 
however, all observables non-invariant under Zjv vanish in the vacuum. 
The action of the symmetry generator does not affect the state |0). Hence, 
the state |S) is locally exactly the same as the state |0) except along the 
boundary C. Therefore, the only regions of space which contribute to 
the overlap are those which lay within one correlation length from the 
boundary. Thus, 

(W(C)) = exp{- 7 P(C)}, (52) 

where P{C) is the perimeter of the curve C. 

The crucial requirement for this argument to hold is the existence of 
a mass gap in the theory. If the theory contains massless excitations the 
factorization of the overlap does not hold, and so in principle even in the 
broken phase the Wilson loop can have a perimeter behavior. This indeed 
is the case in the Abelian theories. 

We now turn to the discussion of low-energy effective theories. This 
will allow us to see an explicit realization of this general argument. 

4.2 The Georgi-Glashow model 

Let us again start with the Georgi-Glashow model. For simplicity all 
explicit calculations in this section will be performed for the SU(2) gauge 
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theory. Generalization to the SU(N) group is not difficult and is discussed 
in Ref. 15. 

Not much has to be done here to parallel the calculation of {V) of 
the previous section. The theory is weakly interacting, and all calcula¬ 
tions are explicit. Choosing the unitary gauge n“ = <5“®, the perturbative 
calculation becomes essentially identical to that in the Coulomb phase of 
QED. The only difference is that the charged matter helds are vectors 
iW^) rather than scalars {(p), but this only enters at the level of the loop 
corrections. The nonperturbative monopole contributions are there, but 
they affect the value of (V) very little, since {V) yf 0 already in perturba¬ 
tion theory. Thus just like in the Coulomb phase of QED, {V) yf 0 and the 
magnetic symmetry is spontaneously broken. The real difference comes 
only when we ask ourselves what is the effective Lagrangian that describes 
the low-energy physics. Here the monopole contributions are crucial, since 
as we have seen before the U(l) magnetic symmetry of QED is explicitly 
(anomalously) broken by these contributions to Zi. Therefore the effec¬ 
tive Lagrangian must have extra terms which reduce the symmetry of the 
Eq. The relevant effective Lagrangian is 

L = d^V*d^V - X{V*V - -b C,{V^ + {y*f) . (53) 

The addition of this extra symmetry breaking term has an immediate 
effect on the mass of the “would be” photon - the phase of V. Expanding 
around {V) = we see that the phase held now has a mass rripf^ = This 
is consistent with the classical analysis by Polyakov Q - the monopole con¬ 
tributions turn the massless photon of QED into a massive (pseudo)scalar 
with an exponentially small mass, niph oc exp{—47rMw/g^}. As a matter 
of fact, for very weak coupling, when the modulus of V can be consid¬ 
ered as frozen, the Lagrangian ( ^^ in terms of the phase x reduces to 
Polyakov’s dual Lagrangian. The exact correspondence between the two 
is discussed in Ref. 15. 

The explicit symmetry breaking causes a dramatic change in the topo¬ 
logically charged (soliton) sector. We know from Polyakov’s analysis that 
the charges in this model are conhned by linear potential, as opposed 
to QED where conhnement is logarithmic. In the effective Lagrangian 
description this is due to the explicit symmetry breaking term. The cru¬ 
cial point is that the vacuum of the theory is not inhnitely degenerate 
{V) = with arbitrary constant x, as in QED, but only doubly degen¬ 
erate {V) = ±p. Thus, the lowest energy state in the nontrivial winding 
sector cannot be a hedgehog. In the hedgehog configuration the field V 
at each point at spatial infinity points in a different direction in the inter¬ 
nal space. This is OK if all these directions are minima of the potential. 
Then the total energy of the configuration comes from the kinetic term, 
and as we saw, is logarithmic. However, now the potential has only two 
minima. Thus, the hedgehog held is far from the vacuum everywhere in 
space. The energy of such a state diverges as the volume of the system, 
E oc g^m^L^. Clearly, to minimize the energy in a state with a nonzero 
winding, the system must be in one of the two vacuum states in as large 
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Figure 2: The string-like configuration of the field V in the state of unit charge 
in the presence of the symmetry breaking terms in the effective Lagrangian. 


a region of space as possible. However, since the field has to wind when 
one goes around the position of a soliton even at an arbitrarily large dis¬ 
tance, V cannot be aligned with the vacuum everywhere at infinity. The 
best bet for a system is, therefore, to choose a string-like configuration. 
Fig. ^ The phase of V(x) deviates from 0 (or vr) only inside a strip of 
width d ~ llrriph stretching from the location of the soliton (charge) to 
infinity. The energy of such a configuration diverges only linearly with the 
dimension L. In fact, a back-of-the-envelope estimate with the effective 
Lagrangian (^^ gives the following energy of such a confining string: 

E oc g^rUphL. 


Clearly the energy of a soliton and an anti-soliton separated by a large 
distance R is E = aR, with the string tension a oc g^niph- 

This is a simple picture of confinement in the effective Lagrangian 
approach in the weakly coupled regime. 

The preceding discussion pertains to confinement of “adjoint” color 
charges. So far we have been considering topological solitons with the 
unit winding, which corresponds to the charge of the massive bosons, 
or “massive gluons.” It should be noted that the notion of the adjoint 
string tension is not an absolute one. So far, our discussion neglected the 
fact that the solitons have a finite core energy and, therefore, in principle 
can be created in pairs from the vacuum. Thus, the soliton-antisoliton 
interaction at a distance R can be screened by creating such a pair, if R 
is big enough. The distance at which the string breaking occurs can be 
estimated from the energy balance between the energy stored in the string 
Es = cTAdjR and the core energy of the soliton-antisoliton pair, which is 
twice the mass of the VF-boson, 2Mw, in our model, 

g^mphR = 2Mw ■ (54) 
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( 55 ) 


The distance at which the string breaks is 


^breaking ^ 


Mw 1 

rriph ' 


Since the width of the string is of the order of Ifniph, and Mw ^ g^ at 
the weak coupling, the length of the string is indeed much greater than 
its width. One can therefore sensibly talk about a well formed adjoint 
confining string. 

As opposed to the adjoint string tension, the concept of the funda¬ 
mental string tension is sharply dehned. This is because the theory does 
not contain particles with fundamental charge and, thus, an external fun¬ 
damental charge can not be screened. To discuss confinement of external 
fundamental charges we have to learn how to deal with half-integer wind¬ 
ings. Imagine adding some extra very heavy fields in the fundamental 
representation to the Georgi-Glashow model. The quanta of these fields 
will carry a half-integer “electric” charge (^) and will be confined with 
a different string tension than To calculate this string tension we 

should consider the Abelian Wilson loop with a half-integer charge. We 
will now do it in the effective theory framework. 

Let us first consider a space-like Wilson loop. As discussed in the 
previous section, this operator is closely related to the generator of the 
magnetic Z 2 . In fact, W(C) is nothing but the operator that performs 
the Z 2 transformation inside the area bounded by the contour C. 

It is straightforward to write down an operator in the effective theory, 
in terms of the field V, that has the same property. 




(56) 


Here S is the surface bounded by the contour C, and P is the momentum 
operator conjugate to the phase of V. In terms of the radius and phase 
of V the path integral representation for calculating the vacuum average 
of this Wilson loop is 

(W(C)) = J Wexpji J d\p^{dp,x- + {d^pf - U{V)^ , (57) 

where U{V) is the Z 2 invariant potential of Eq. (^). The external current 
(x) does not vanish only at points x which belong to the surface S and 
is proportional to the unit normal to the surface S. Its magnitude is 
such that when integrated in the direction of n it is equal to tt. These 
properties are conveniently encoded in the following expression: 


dxpj^ix) = ■ 


i{T,C). 


(58) 


Here T is an arbitrary closed contour, and n(T, C) is the linking number 
between two closed curves T and C. 


^Note that here we are dealing with the path integral representation, and thus the contour 
C and the surface S are embedded into a three dimensional Euclidean space. The linking 
number between two curves is also defined in three dimensions. 
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Figure 3: The external current which creates a pair of static fundamental 
charges in the effective Lagrangian description. 


This path integral representation follows immediately if we note that 
the conjugate momentum P is P{x) = 2p^{x)doxix) some fixed time 
t. This accounts for the linear in d^x term in the exponential (p^. The 
constant term arises due to the standard integration over the conjugate 
momenta in passing to the path integral representation. 

The path integral representation was constructed for the spatial Wil¬ 
son loops. However, the expression (|^) is completely covariant, and in 
this form it is valid for time-like Wilson loops as well. It is important to 
note that, although the expression for the current depends on the surface 
S, the Wilson loop operator depends only on the contour C that bounds 
this surface. A simple way to see this is to observe that a change of vari¬ 
ables X ^ X + in the volume bounded hy S + S' leads to the change 
in Eq. The potential is not affected by this change since it 

is globally Z 2 invariant. Therefore, the operators defined with S and S' 
are completely equivalent. 

To calculate the energy of a pair of static fundamental charges at the 
points A and B we have to consider a time-like fundamental Wilson loop 
of an infinite time dimension. This corresponds to time-independent jfj, 
which does not vanish only along a spatial curve G (in the equal-time 
cross section) connecting the two points and pointing in the direction 
normal to this curve, FigM The shape of the curve itself does not matter, 
since changing the curve without changing its endpoints is equivalent to 
changing the surface S in Eq. 

In the classical approximation the path integral ([stI) is dominated by a 
static configuration of V. To determine it we have to minimize the energy 
on the static configurations in the presence of the external current j^. 
The qualitative features of the minimal energy solution are quite clear. 
The effect of the external current is to flip the phase of 1/ by tt across 
the curve G, as is expressed in Eq. (^). Any configuration that does not 
have this behavior will have the energy proportional to the length of G 
and to the ultraviolet cutoff scale. Recall that the vacuum in the theory at 
hand is doubly degenerate. The sign change of V transforms one vacuum 
configuration into another. Therefore, the presence of requires that 
on the opposite sides of the curve G, immediately adjacent to G, there 
should be different vacuum states. It is clear, however, that far away from 
G in either direction the held should approach the same vacuum state, 
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Figure 4: The minimal energy configuration of V in the presence of a pair of 
fundamental charges. 



otherwise the energy of the configuration diverges linearly in the infrared. 
The phase of V has to make half a wind somewhere in space to return to 
the same vacuum state far below G, as the vacuum state that exists far 
above G. If the distance between A and B is much larger than the mass 
of the lightest particle in the theory, this is achieved by having a segment 
of a domain wall between the two vacua connecting the points A and B. 
Clearly to minimize the energy, the domain wall must connect A and B 
along a straight line. The energy of such domain wall is proportional to 
its length, and, therefore, the Wilson loop has the area law behavior. The 
minimal energy solution is schematically depicted in Fig. 

We see that the string tension for the fundamental string is equal 
to the tension of the domain wall which separates the two vacua in the 
theory. This relation has been discussed a long time ago Q by ’t Hooft. 
Parametrically, this string tension is clearly the same as the adjoint one, 

cr/ oc g^ruph , (59) 

although the proportionality constant is different. We will briefly discuss 
the relation between the adjoint and the fundamental string tensions in 
the next subsection. 

Note that the fundamental string is an absolutely stable topological 
object in the Z 2 invariant theory: the domain wall. It can not break, if 
one makes the distance between the two charges larger. In the effective 
theory it is also obvious since there are no point-like (particle-like) objects 
in the theory on which the domain wall can terminate since there are no 
dynamical objects with half-integer winding number. 
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We interpreted this calculation as the calculation of the potential be¬ 
tween the fundamental adjoint charges - the time-like Wilson loop. How¬ 
ever, in the Euclidean formulation there is no difference between the time¬ 
like and space-like Wilson loops. Interpreted in this way the calculation 
becomes a technical illustration to the argument given in the previous 
subsection: space-like Wilson loops have the area law if the magnetic Z 2 
symmetry is spontaneously broken. 

4.3 Gluodynamics 

In the weakly interacting case the effective low-energy Lagrangian can 
be derived as explained above in perturbation theory, combined with the 
dilute monopole gas approximation. The more interesting regime is cer¬ 
tainly that of strong coupling, which is essentially pure Yang-Mills theory. 
The luxury one has in 2-1-1 dimensions is that the weak and strong cou¬ 
pling regimes are not separated by a phase transition.This means that 
whatever global symmetries the theory has, their realization must be the 
same in the weakly coupled and strongly coupled vacua. The existence of 
the Zjv symmetry is an exact statement not related to the weak coupling 
limit. Therefore, it is natural to expect that this symmetry must be non- 
trivially represented in the effective low-energy Lagrangian. It is plausible 
then that the low-energy dynamics at strong coupling is described by the 
same effective Lagrangian which encodes the spontaneous magnetic sym¬ 
metry breaking. Certainly, the values of the coupling constants will be 
different in two regimes, but the qualitative behavior should be similar. 

Strictly speaking, this is an assumption rather than theorem. That is 
where the difference between continuous and discrete symmetries comes 
in. Were magnetic symmetry continuous (like in the Abelian case), its 
spontaneous breaking would unambiguously determine the structure of 
the effective Lagrangian, whether the theory was weakly or strongly cou¬ 
pled. With discrete symmetries this is not necessarily the case. It could 
happen that even though the symmetry is broken, the “pseudo-Goldstone” 
particle is so heavy that it decouples from the low energy dynamics. For 
this to happen, though, the symmetry breaking would have to occur at 
a very high energy scale. In gluodynamics this is very unlikely, since the 
theory has only one dynamical scale. In fact, as was discussed in the 
beginning of this section, the fundamental string tension determines the 
scale at which Zn is broken. Wilson loops of linear size Z < do 

not distinguish between confining and nonconfining behavior, and hence 
between the broken and the unbroken Zjv. The scale of Zjv breaking is 
which is precisely a natural dynamical scale of QCD. 

Therefore, generally speaking, we expect that the pseudo-Goldstone 
meson stays among the low-energy excitations not only in the weakly 
coupled limit but also in pure gluodynamics. If this is the case the degrees 
of freedom that enter the effective Lagrangian in the weakly coupled phase 
also interpolate actual low-energy physical states of the strong coupling 
regime. That is to say, the radial and phase components of the vortex 
field V must correspond to the lightest glueballs of pure SU(Y) Yang 
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Mills theory. We can check whether this is the case by considering the 
lattice gauge theory data |^] on the spectrum. The radial part of V is 
obviously a scalar and has quantum numbers 0’*'^. The quantum numbers 
of the phase are easily determined from the definition ( [l4[ ) . Those are 0 . 

The spectrum of pure SU(A^) Yang Mills theory in 2+1 dimensions was 
extensively studied recently ||m on the lattice. The two lightest glueballs 
for any N are found to have exactly these quantum numbers. The lightest 
excitation is scalar, while the next one is a charge-conjugation-odd pseudo¬ 
scalar meson with the ratio of the masses roughly u mpjms = 1.5 for any 
N. 

The situation, therefore, is likely the following. The low-energy physics 
of the SU (2) gauge theory is always described by the effective Lagran gian 
d^). In the weak coupling regime the parameters are given in Eqs. 
and Here the pseudo-scalar particle is the lightest state in the spec¬ 
trum while the scalar is the first excitation. The pseudo-scalar meson 
presents an “almost massless photon.” The scalar meson is the massive 
Higgs particle. Moving towards the strong coupling regime (decreasing 
the Higgs VEV) one increases the pseudo-scalar mass while reducing the 
scalar mass; the parameters of the effective Lagrangian change accord¬ 
ingly. The crossover between the weak and strong coupling regimes occurs 
roughly where the scalar and the pseudo-scalar become degenerate. At 
strong coupling the degrees of freedom in the effective Lagrangian are the 
two lightest glueballs. They are still collected in one complex held which 
nontrivially represents the exact Z 2 symmetry of the theory (or Zn for 
SU(N)). 

Of course the spectrum of pure Yang-Mills theory, apart from the 
scalar and the pseudo-scalar glueballs, contains many other massive glue- 
ball states and those are not separated by a large gap from the two lowest 
ones. Application of the effective Lagrangian in the strong coupling regime 
has to be considered as qualitative rather than quantitative. At this qual¬ 
itative level though, linear conhnement is an immediate property of the 
effective Lagrangians of this type. 

The fact that the masses of the scalar and pseudo-scalar particles 
are interchanged in gluodynamics relative to the weakly coupled regime 
leads 0 to some interesting qualitative differences. In particular, the 
structure of the confining string and the interaction between the strings 
differ in some important ways. Let me briefly discuss this aspect. 

In the weakly coupled regime the phase of V is much lighter than 
the radial part. A cartoon of the fundamental string in this situation is 
depicted in Fig. ^ 

The radial part p being very heavy practically does not change inside 
the string. The value of p in the middle of the string can be estimated 


^Actually this pattern is firmly established only for N > 2. At N = 2 the mass of the 
pseudo-scalar meson was not calculated. ]l0[ The reason is that it is not clear how to construct 
a charge-conjugation-odd operator in pure gauge SU(2) lattice theory. So, it is possible that 
the situation at A = 2 is non-generic in this respect. In this case our strong coupling picture 
should apply at A > 2. 
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Figure 5: The structure of the string (domain wall) in the regime when the 
pseudo scalar is lighter than scalar, m < M. 


from the following simple argument. The width of the region where p 
varies from its vacuum value p to the value po in the middle, is of the 
order of the inverse mass of p. In terms of the energy per unit length this 
variation costs 

2 

apr^ M{p-pof+ x'^pl, (60) 

where a; is a dimensionless number of order unity. The first term is the 
contribution of the kinetic term of p while the second contribution comes 
from the interaction term between p and x due to the fact that the value of 
y in the middle of the string differs from its vacuum value. Our notations 
are such that m is the mass of the pseudo-scalar particle and M is the 
mass of the scalar one. Minimizing this with respect to po we find 

po = (61) 

Thus, even in the middle of the string the difference in the value of p and 
its VEV is of the second order in the small ratio m/M. Correspondingly, 
the contribution of the energy density of p to the total energy density is 
also very small, 

~ • (62) 

This is to be compared with the total tension of the string which is con¬ 
tributed mainly by the pseudo-scalar phase y, 

~ mp^ . (63) 

We obtain this by estimating tie kinetic energy of x on a configuration 
of width 1/m where x changes u by an amount of order 1. 

^The fact that the heavy radial field p practically does not contribute to the string tension 
is natural from the point of view of decoupling. In the limit of infinite mass p should decouple 
from the theory without changing its physical properties. It is very different,however, from 
the situation in superconductors. In the superconductor of the second kind, where the order 
parameter field is much heavier than the photon ( k > -^) the magnetic field and the order 
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Figure 6: The structure of the string (domain wall) in the regime when the 
pseudo-scalar is heavier than the scalar, m > M. 


Let us now consider the domain wall (or fundamental string) in the 
opposite regime, when the mass of the scalar meson is much smaller than 
that of the pseudo-scalar one. Now, the profile of the fields in the wall is 
very different. A cartoon of this situation is given in Fig. ^ 

We will use the same notations, denoting the mass of the pseudo-scalar 
by m and the mass of the scalar by M, but now m ';$> M. Let us again 
estimate the string tension and the contributions of the scalar and of a 
pseudo-scalar to it. The width of the region in which the variation of p 
takes place is of the order of its inverse mass. An estimate of the energy 
density of the p field is given by the contribution of the kinetic term 

(jp ~ M(^ - po)^ . (65) 

The width of the region in which the phase x varies is ~ 1/m. In this 
narrow strip the radial held p is practically constant and is equal to po. 
Hence, the kinetic energy of x contributes 

~ mpo . (66) 


Minimizing the sum of the above two contributions with respect to po we 
hnd 


Po 


M 


— p <C p . 
m 


(67) 


parameter give contributions of the same order (up to logarithmic corrections 0(log k)) to the 
energy of the Abrikosov vortex. This is the consequence of the fact that the order parameter 
itself is forced to vanish in the core of the vortex, and therefore even though it is heavy, its 
variation inside the vortex is large. An even more spectacular situation arises if we consider 
a domain wall between two vacuum states in which the heavy field has different values. ]l l| In 
this situation the contribution of the heavy field (j) to the tension would be 

f^heavy = M (64) 

where Acf) is the difference in the values of (p on both sides of the wall. For a fixed the 
energy density diverges when (p becomes heavy. In the present case this does not happen since 
the two vacua which are separated by the domain wall differ only in VEV of the light field x 
and not the heavy field p. 
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Also 


Mo o 

a^c ~ —M/i , ap = Op M ^ . (68) 

m 

Now the radial field is very small in the core of the string. The energy 
density is contribuit|ed almost entirely by the scalar rather than by the 
pseudo-scalar field.l3 

The extreme situation m ^ M is not realized in the non-Abelian 
gauge theory. From the lattice simulations we know that in actuality 
even in the case of pure Yang-Mills the ratio between the pseudo-scalar 
and scalar masses is about 1.5 - not a very large number. The analysis 
of the previous paragraph does not reflect the situation in the strongly 
coupled regime of the theory. Rather we expect that the actual profile of 
the string is somewhere in between Fig. ^ and Fig. ^ although somewhat 
closer to Fig. The widths of the string in terms of the scalar and 
pseudo-scalar fields are of the same order, although the scalar component 
is somewhat wider. The same is applicable to the contribution to the 
string tension. Both dueballs contribute, with the scalar contribution 
being somewhat larger .□ 

The interaction between two domain walls (confining fundamental 
strings) in the two extreme regimes is also quite different. In the weakly 
coupled region we can disregard the variation of p. For two widely sepa¬ 
rated parallel strings the interaction energy comes from the kinetic term 
of X- This obviously leads to a repulsion, since for both strings in the 
interaction region the derivative of the phase is positive. On the other 
hand, if the pseudo-scalar is very heavy the main interaction at large sep¬ 
aration is through the “exchange” of the scalar. This interaction is clearly 
attractive, since if the strings overlap, the region of the space where p is 
different from its value in the vacuum is reduced relative to the situation 
when the strings are far apart. 

Therefore, the situation is very similar to that in superconductivity. 
The confining strings in the weakly coupled and strongly coupled regimes 
behave like the Abrikosov vortices in the superconductor of the second and 
first kind, respectively. This observation has an immediate implication 
for the string tension of the adjoint string. As we have discussed in the 
previous section the phase x changes from 0 to 27r inside the adjoint string. 
The adjoint string can be pictured as two fundamental strings running 
parallel to each other. In the weak coupling regime the two fundamental 
strings repel each other. Therefore, two fundamental strings within the 
adjoint string will not overlap, and the energy of the adjoint string is twice 
the energy of the fundamental one, 


cTAdi = 2crF . (69) 

®Again, this is in agreement with decoupling. The heavier field does not contribute to the 
energy, even though its values on the opposite sides of the wall differ by 0(1). Its contribution 
to the energy is suppressed by the factor Pq which is very small inside the wall. 

®One also expects similar non-negligible contributions from higher mass glueballs which 
are not taken into account in our effective Lagrangian framework. 
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In the strongly coupled region the situation is quite different. The 
strings attract. It is clear that the contribution of p to the energy will be 
minimized if they overlap completely. In this case the contribution of p in 
the fundamental and adjoint strings will be roughly the same. There will 
still be a repulsion between the pseudo-scalar cores of the two fundamental 
strings, so presumably the core energy will be doubled. In this way we 
arrive at an estimate 


aAd^ ^ (IF + 0{ —a f) ■ (70) 

Again the situation is more complicated in the limit of pure Yang- 
Mills. The scalar is lighter, therefore the interaction at large distances is 
attractive. However, the pseudo-scalar core size and its contribution to 
the tension are not small. In other words. M/m in Eq. is a number 
of order one. 

Thus, in broad terms the relation between the weakly and strongly 
coupled confining theories is similar to the relation of the superconductors 
of the second and first kind. 


5 Fundamental Matter. 

The discussion in the preceding sections pertained only to theories with¬ 
out matter fields in fundamental representation. Obviously it is very 
important to understand what becomes of the magnetic symmetry and 
the associated order parameter if fundamental matter is present. In this 
section we address this questionjl^. 

For simplicity we will consider scalar quarks. The Lagrangian of the 
theory we are interested in is 

L = -iTrE^ + |L>^<f>|^-M^<l>*$ (71) 

where the scalar field "h™ transforms according to the fundamental repre¬ 
sentation of the SU{N) color group. 

5.1 Zn as a local symmetry. 

First thing to note is that the fundamental Wilson loop still commutes 
with the Hamiltonian. This is obvious, since the extra terms in the Hamil¬ 
tonian in the presence of the quarks do not involve electric field operator, 
but only vector potential. The Wilson loop commutes with the vector 
potential, and thus with the additional terms in the Hamiltonian. The 
theory thus still has the magnetic Zn symmetry. 

This may seem somewhat surprising at first sight. We have shown 
above that spontaneous breaking of the magnetic Zn implies the area 
law for the Wilson loop, and conversely the perimeter law of W implies 
unbroken Zn- In the theory with fundamental charges the Wilson loop 
is known to have perimeter law due to breaking of the confining string 
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at any finite value of the fundamental mass M. We might then conclude 
that the magnetic Zn is restored at any, arbitrarily large value of M but 
is broken at M —> oo. The common lore is that the Zn breaking phase 
transition for > 2 is first order, and this then implies a discontinuous 
behaviour of the theory in the infinite mass limit. This of course is com¬ 
pletely counterintuitive and in fact plain wrong. The caveat in this line 
of reasoning is the following. The relation between the behaviour of the 
Wilson loop and the mode of the realization of the magnetic symmetry 
hinges crucially on the existence of a local order parameter of the mag¬ 
netic Zn . In the absence of such an order parameter it is not true that the 
Wilson loop locally changes the quantum state inside the loop and this 
invalidates the whole argument. In particular in the absence of a local 
order parameter, the Zn symmetry can be spontaneously broken but the 
Wilson loop can have a perimeter law. 

In fact it is easy to see that for any hnite M the magnetic Zn does 
not have a local order parameter. The only candidate for such an order 
parameter is the vortex operator V{x) defined in eq.(^^, since it has to 
be local also relative to the purely gluonic operators. However in the 
presence of fundamental quarks the operator V(x) is not local anymore. 
To see this consider the dependence of Vc'(*) on the curve C which enters 
its dehnition. As before the operators Vc and Vc' are related by 

Vcix) = Vc'{x)exp{^ j d^xTrdiYE'} (72) 

where S is the area bounded by C — C'. As before, due to the Gauss’ law 
the integral in the exponential is equal to the total hypercharge in the area 
S. However the hypercharge of fundamental quarks has eigenvalues Eg/2. 
The extra phase factor is therefore not unity anymore but can rather take 
values exp{27ri/A’} depending on the state and the choice of the contour 
C. 

The status of the magnetic Zn is thus quite different in a theory with 
fundamental quarks - it does not have a local order parameter. Neverthe¬ 
less it is clear that at least as long as the mass M is large M/g^ > 1, the 
relevant degrees of freedom for the effective infrared dynamics should still 
be the vortices V and the main factor which determines their dynamics 
should still be the magnetic Zn- At large M the dependence of V on the 
curve C is weak, since the probability of the vacuum fluctuations which 
involve fundamental charges is small. The probability of appearance of a 
virtual qq pair separated by a distance I is suppressed by the exponential 
factor exp{—MZ}. The typical distance scale for the ’’glueball” physics 
is 1/g^. Thus at these distances such fluctuations are unimportant and 
should not affect much the dynamics. Things are different however if one 
is interested also in the baryonic sector of the theory. The baryons are 
necessarily heavy and in order to be able to discuss their structure we 
must understand the main dynamical effects also at shorter distances. 

Therefore our aim now is to understand what is the main effect of the 
non locality discussed above on the dynamics of magnetic vortices. 
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The situation we have just described - a symmetry without a local 
order parameter - is not exceptional in quantum field theory. This is pre¬ 
cisely the property of the global part of any Abelian gauge group. Con¬ 
sider for example quantum electrodynamics. The global electric charge is 
of course a physical gauge invariant charge with the corresponding gauge 
invariant local charge density. 


0 = 



(73) 


Nevertheless there is no local operator that carries this charge. This is a 
direct consequence of the Gauss’ law 


diEi = gp (74) 

Any physical, gauge invariant operator that carries Q must also carry the 
long range electric field, which can not fall off faster than a power of the 
distance. The gauge invariant QED Lagrangian is written in terms of 
’’local” charged fields (f). But appearances are deceptive: these fields are 
not gauge invariant, and therefore not physical. A gauge invariant charged 
field can be constructed from 4> by multiplying it by a phase factor 

<t>phys{x) = (t>exp{ig j d^yei{x - y)Ai{y)} (75) 

with the c-number field a satisfying 

diei = 5^{x) (76) 

For any a satisfying this condition the operator 4>vhys is gauge invari¬ 

ant, and therefore physical. It is however necessarily nonlocal. Different 
choices of Ci define different gauge invariant operators and correspond to 
different gauge fixings. Thus for ei{x) = Xijx^ the field (j)phya is the field 
cj) in the Coulomb gauge, while ei{x) = SiiS(x 2 )ff(xi) corresponds to the 
axial gauge Ai = 0, and so on. Different definitions of (jiphya differ from 
each other by a phase factor, which is precisely the gauge ambiguity of 
the original field (f>. 

The U{1) gauge group is the most natural Abelian gauge symmetry to 
consider in continuum field theory. One can however also consider discrete 
groups like Zjvj^. In this case again no local operator that carries the 
global Zn charge exists. The various gauge invariant charged operators 
are nonlocal and differ from each other by a local Zn - valued phase. 
These different operators again correspond to different gauge fixings of 
the local Zn group. 

This is precisely the structure that emerged in our discussion in the 
earlier part of this section. We have 

1. Global magnetic Zn symmetry generated by the fundamental Wil¬ 
son loop. 

2. The set of nonlocal vortex operators Vc{x), which all carry the Zn 
charge and differ from each other by a Zn valued phase factors. 
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It is very suggestive therefore to think about Vc(x) as of different 
gauge fixed versions of a field charged under local Zm- This leads us to 
expect that the low energy theory we are after should be a Zn gauge 
theory of the magnetic vortex field V. 

In fact coming back to the discussion in the beginning of this section, 
we see that from this vantage point it is obvious why the Wilson loop 
has a perimeter law, even if the global Zjv is broken spontaneonsly. The 
action of a Wilson loop of a finite size inside the contonr is a gauge trans¬ 
formation. Thus in physical terms locally inside the contonr the new state 
is the same as the old one and so the overlap between the two locally is 
unity. The only nontrivial contributions to the overlap come from the 
region close to the contour, thus giving the perimeter law. 

We will now construct the effective low energy locally Zm invariant the¬ 
ory, and discuss its relation to the original QCD Lagrangian. Although it 
is possible to give a more ‘analytic’ derivation of the effective Lagrangian, 
we prefer to write it down directly guided by the previous discussion. We 
will then explain the physical meaning of the various fields that appear in 
it. 


5.2 Gauging the Wilson loop. 

The easiest way to construct a Zn gauge theory in the continuum is to 
consider a C/(l) gauge theory with the Higgs field of charge N which has a 
large expectation value |^. Consider therefore the following Lagrangian: 

L = + + \{d^-ib^)U\^ 

- \{V*V - i/f+ V*^U)-X{U*U-uf . (77) 

Here = d^bv — dvbfj,. We take the parameters such that, A >> A and 

» fj.^. 

The ‘Higgs’ field U has a large expectation value and breaks the U{1) 
gauge symmetry down to its Zn subgroup V{x) —> }V{x). 

Below the scale determined by the expectation value u, the field U is 
practically frozen and its fluctuations are unimportant. In this regime the 
model indeed describes the locally Zn invariant theory. The global part 
of the gauge group is our Zn magnetic symmetry generated by the Wilson 
loop. The larger 17(1) gauge structure at this point is just an auxiliary 
trick which enables us to write down a discrete gauge theory in continuous 
notations. We will however see later that it does in fact has a real physical 
meaning of its own and arises naturally in the effective theory. 

Let us first discuss how this Lagrangian reduces to the effective La¬ 
grangian of the pure Yang-Mills theory eq.jjs^ in the limit of the large 
quark mass. In this limit not only the field U must decouple, but also the 
gauge interactions of the field V must vanish. There is another consistency 
requirement. In the limit of zero gauge coupling —> 0 the gauge 17(1) 
symmetry becomes global 17(1) and is broken due to non vanishing expec¬ 
tation value of U. The spectrum therefore contains a massless Goldstone 
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boson. This Goldstone boson is of course the longitudinal component of 
the gauge field b^. There is however no such massless particle in the pure 
gluodynamics nor in the effective Lagrangian eq.(^). This means that 
the couplings in eq.(|^) should depend on the quark mass in such a way 
that the vector particles remains heavy for any finite M and its mass goes 
to zero very sharply only in the limit when it is completely decoupled. In 
terms of the Goldstone boson couplings it means that must be larger 
than any scale relevant to the dynamics of the vortex field V. All these 
conditions can be met by choosing for example 


2 

e = 


M’ 


u = xM, 




(78) 


With this choice the mass of the vector boson stays finite 

as M —> oo and can be arbitrarily large. The Goldstone boson in the 
decoupling limit has an infinite f-n and is completely decoupled just as the 
’’invisible axion”. 

How do the basic fields present in the effective Lagrangian arise in 
the fundamental theory eq.(^^? To understand this, let us start with 
considering the symmetries of the theory. 

The Lagrangian eq.([7^ has two global U{1) symmetries. One is the 
global part of the local U{1) with the conserved current ej J = dvfu/j ■ The 
other current, conserved by virtue of the homogeneous Maxwell equation 
is /m = ^iivxfvX- The two charges have quite different nature. The second 
one, the dual magnetic flux 


<E>£) 



(79) 


has a local order parameter. It can be constructed in a way similar to the 
vortex field in QED §.§■ The global 17(1) gauge charge, on the other 
hand does not have a local order parameter, as discussed earlier. 

Both these conserved currents also exist in the fundamental theory. It 
is fairly straightforward to identify them. The QGD Lagrangian has one 
obvious global 17(1) charge - the baryon number. This charge has local 
order parameters - gauge invariant baryon fields of QCD, and is therefore 
identified with the dual magnetic flux 

Qb = <i>D (80) 

The second charge can be expressed in terms of the spatial current 
components of the first one 


Q 


T 



xdi 


r fj 


(81) 


It is thus the vorticity associated with the baryon number current. 

The vortex operator V which appears in eq. ([tt]) is not gauge invariant 
and is only physical after complete gauge fixing of the 17(1) gauge group. 
After such a gauge fixing (which amounts to multiplying V by an operator 
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valued phase) the Gauss law requires that on the physical states the phys- 
ical operator V carry the charge Qt- Due to the identification eqs.( p^ , pi[ ) 
this leads to an interesting conclusion that any physical operator V in the 
effective theory creates a vortex of the original baryon number current. 

This somewhat unexpected conclusion is in fact quite natural for an 
eigenoperator of the magnetic Zjv. Consider the vortex operator V defined 
by eq. (p^ . In QCD just like in pure gluodynamics, it has an alternative 
representation of an operator of the singular gauge transformation of the 
form eq.(Q 

Vc(a;)=exp{^ J (YE j(y)) + 9(x - y)Jo (y)} ( 82 ) 

The only difference is that now is the hypercharge operator due to 
both, gluons and fundamental quarks 


Jo"" = ig \TrY[Ai, Ei] + - -h^n;) 


(83) 


Consider the action of this operator on a quark field "I?. The transformed 
quark field >!>' = is a gauge transform of $ everywhere except 

along the branch cut of the function 9. Across this cut the phase of 
is discontinuous - it jumps by 2-n/N for all color components of Due 
to this discontinuity the baryon current - the global U{1) current of "I? - 
does not vanish at points along the cut. Calculating explicitly the action 
of V on the baryon number current we find 


V^jf{x)Vc = = jf {x)+‘^nt {x)5{x G C)-I>*cE-(a:) 

(84) 

where nf (x) is a unit vector normal to the branch cut C at the point x. 

It is natural to define the local vorticity associated with the baryon 
number as 


Pt = mjdi[ 




(85) 


The vortex operator V therefore creates a vortex of baryon number current 
with fractional vorticity 2n/N 


VMx)pT{y)Vc{x) = priy) + - y) (86) 


The operator U due to the Gauss’ law also carries baryon vorticity. 
Since its gauge coupling is N times the coupling of V, it creates one unit 
of vorticity. 

This simple exercise also helps us to identify the value of the gauge cou¬ 
pling constant in the effective theory. Comparing eq.(^^ with eq.(^l|) 
we find 

oc <&*<& (87) 

The same relation is obtained by comparing the current algebra in the 
fundamental and the effective theories. The commutator of the baryon 
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charge density with the baryon current density in the fundamental theory 
is 

[Jo{x),jf{y)] = i$*^diS^{x-y) (88) 

In the effective theory using the canonical commutators that follow from 
eq. ( 0 ) and the identification eq.(^^ we find 

2 

[Joix),jf{y)]=i^diS'^{x-y) (89) 

Again we deduce eq. (^7|) . The operator on the right hand side of this 
equation in the effective theory is indeed a constant. The effective theory 
should be valid at long distances. In this regime in the leading order in 
the derivative expansion the operator $*<!> should be approximated by its 
expectation value. Taking into account fluctuations of is tantamount 
to including higher derivative terms in the effective Lagrangian eq. ( 0 )- 
At this, higher order in derivative expansion the gauge coupling constant 
would become a dynamical field. While this is perfectly legitimate, it is 
beyond our present framework. 

When calculating the expectation value of >!>*$ we should remember 
that it has to be calculated with the cutoff A appropriate for the effective 
theory. This cutoff must be above the characteristic scale of the pure 
gluodynamics (determined by the string tension) but below the heavy 
quark mass. With this in mind we get 


< <&*'!' >A= 


d^p 


8n^ (p2 + M2)i/2 


M 


(90) 


So that finally 


2 


e 



(91) 


which is consistent with the expected scaling in the large mass limit 

eq.(0. 


5.3 The baryon and the bag. 

Having understood the origin of the fields and the symmetries in the 
effective theory, we would like to see how it encodes the qualitative features 
of the low energy QCD physics. Since we are considering a heavy quark 
theory, below the fundamental mass scale the spectrum should be the 
same as in pure gluodynamics. We have already seen that, in the infinite 
mass limit, the effective theory reduces to that of (^^. Indeed, even at 
finite but large M, this is the case at low energies. Since the VEV of the 
field U is large, we can impose the unitary gauge condition on it. In this 
unitary gauge the phase of U disappears. The modulus of U is very heavy, 
and so is the vector field 6^. Thus at low energies we recover the effective 
theory of pure glue sector. There is however a set of configurations, on 
which the unitary gauge can not be imposed. Those are configurations in 
which U vanishes at some points in space. Indeed it is these configurations 
that are important for the baryonic sector. Recall that 17 is a vortex of 
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baryon number current. Thns, by dnality one expects that the baryon 
charge is associated with the vortex configuration of the field U. In the 
core of the vortex the field U has to vanish, and so the nnitary gange is 
not admissible. Thus in the baryon sector we can not think of U as frozen 
at its expectation value and instead have to treat it as a dynamical field. 

That the baryon does indeed carry vorticity of the field U, can be seen 
by the following simple argument. The baryon nnmber is represented in 
the effective theory by the dual magnetic flux. The baryonic state must 
therefore be the dual magnetic vortex. Such a vortex in a nonsingular 
gauge has a vector potential of the form bt — To have a finite 

energy it must be accompanied by the winding of the phases of both V 
and U. Since U carries N times the charge of V, the only states that are 
allowed energetically are those that carry N vortices of U. 

This is natural in view of the ‘dual’ relation between the effective 
and the fundamental theories. The field U is the ‘vortex’ dual to the 
fundamental quark. We thus expect that its elementary vortex would 
represent the fundamental quark itself, and so finite energy states must 
contain N such elementary vortices. A single vortex must be confined. 
In the same sense, V is ‘dual’ to the adjoint gluon field. The elementary 
vortex of V is then, in a sense, the ‘constituent’ gluon. In fact, such single 
gluon should not exist as a finite energy state either, and we expect it also 
to be confined. The single vortex of V should therefore bind either with 
the anti-vortex of the same type, or with N vortices of U. The former 
type of state is a glueball, and exists in the pure glue theory , while the 
latter type is a baryon. 

Interestingly enough, this line of reasoning leads us to expect that the 
baryon must have a bag-like structure. Namely, the quarks are bound to 
the V field vortex. Inside this vortex, the value of V is small - it vanishes 
in the middle, and than rises quite slowly (relative to the scale of 1/M) 
towards the edges. Recall that V — 0 corresponds to a non confining 
state The quarks are therefore sitting in the ‘perturbative’ region of 
space - where there are no confining forces. Only when they separate 
far from each other - into the region with non vanishing V, the linear 
potential pulls them inside again. 

Let us look at this more carefully. Consider for simplicity the Z 2 
symmetric case N = 2. The baryon is the dual vortex with dual magnetic 
flux 2-k je. Far from the vortex core, the field configuration is pure gauge, 
with the phases of V and U following the vector potential: 

h = V{x) = U{x) = . (92) 

x-^ 

The parameters of the model are such that the field V is much lighter 
than both U and bi. Thus the size of the vortex core of V is large - of the 
order of the inverse glueball mass. The two [/-vortices which have a very 
small size core, sit inside this core. Since the length associated with the 
dual magnetic field is much smaller than the core size of the V - vortex, 
the dual flux is concentrated on the U - vortices. From the low energy 
point of view, the picture is that two point-like magnetic vortices sit inside 
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Figure 7: The ‘baryon’ configuration for Z 2 . Arrows represent the direction 
(phase) of V, and the big dots at A and B correspond to the positions of the U 
field vortices. 


a soft core of a F field vortex. The field configuration looks roughly as 
depicted on figure 7. 

The magnetic flux is concentrated at the points A and B. The phase 
of the field [/ follows the variation of the vector potential very closely. 
The phase of V is also trying to do that, but it can not quite follow it 
all the way, since on the line between the two vortices it would have to 
be discontinuous. The most important energy contribution (apart form 
the core energy of small vortices) comes therefore from the vicinity of this 
line. The phase of V obviously has to interpolate across this line between 
the values 0 and tt. Since the modulus of V is not extremely rigid, it will 
be smaller along this line than in the immediate neighborhood. Both, the 
variation of the modulus and the phase of V along the line connecting 
the small vortices contribute to the energy which is clearly linear in the 
distance \ A — B\. 

To study the structure of the baryon in the effective theory more quan¬ 
titatively let us fix the gauge such that 

U{x) = Mexp[i0(a;)] 

6{x) = arctan(— - —)-|-arctan(— - —). (93) 

X Q> X — 1 “ CL 

This is a valid gauge choice for the configuration with two vortices of unit 
winding at the positions (a,0) and (—a,0). 
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Due to the conditions on the parameters of our model, the vector field 
6^ follows the phase of U in the whole space. Thus 


bo = 0, bj 


D'fc 


(a; — 

\x — a;d)| 


Ife 

T 


+ D'fc 


(x - 
\x — 


(94) 


where: xi = x, X 2 = y, and I, r are the positions of the left and right 
vortices, respectively: = (—a,0), = (a,0). 

There are two interesting limiting situations. The first is when the 
distance a is larger than the dynamical distance scales and 

which determine the masses of the glueballs. The second is the 
reverse, that is when the two vortices are sitting well inside the glueball 
correlation length. Let us look at them in turn. 

When the distance between the vortices is large we expect the potential 
between them to be linear with the string tension calculated in the theory 
without the dynamical U field. We will study the interquark potential in 
our low energy theory in the semiclassical approximation. 

To find the minimal energy configuration, we have to solve the classical 
equations of motion for the field D at a fixed configuration of U and 6; 
given by eqs.( |93| , |9^ ). Let us concentrate on the points which are close to 
the X - axis, with |x| << a. The main contribution to the energy comes 
from this region of space. In this region the vector potential bi vanishes, 
and the phase of the field U is zero. Thus the equations of motion for the 
field V are the same as in the pure gluodynamics. Also as long as we stick 
to this region, the configuration of V is translationally invariant in the x 
direction. What determines the energy then are the boundary conditions 
on the field V. In this configuration clearly the phase of the field V is n/2 
far above the x axis and —7r/2 far below the axis, fig.7. Thus both, the 
equations and the boundary conditions are precisely the same as for the 
domain wall separating the two degenerate vacua in the effective theory 
of pure gluodynamics (e^ =0, A —> oo). There is an extra contribution 
to the energy that comes from the region of space close to the points A 
and B. But this energy does not depend on a and is subleading for large 
a. The rest of the space does not contribute to the energy, since the field 
configuration there is pure gauge. 

Thus as expected the energy in this regime is i? = aa, where a is 
the domain wall tension (fundamental string tension) calculated in pure 
gluodynamics 

It is also interesting to consider the opposite situation, when the dis¬ 
tance between the quarks is smaller than the glueball correlation length. 
This is the regime in which we do not expect to see any stringy structure. 
Instead we can ask whether the lowest energy configuration has any re¬ 
semblance to a bag. To study this question we take the limit a ^ 0. In 
this case the phase of the field V will follow the phase of U in the whole 
space 



u* ■ 


(95) 


Given this condition, only the variation of the radial component of V 
has to be determined. Since the problem has rotational symmetry, the 
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equation of motion for the modulus p becomes 


d?p 

dr'^ 


— 2rip + 2Xp = 0 


(96) 


with p = — ^u, which is assumed to be positive throughout this 

paper. The relevant boundary condition in this case is that at infinity 
p approaches its vacuum value p ,—.oo ^ v while in the vortex core it 
vanishes p(0) = 0. With these boundary conditions, eq. (^) has the 
familiar form of the static kink equation in 1 + 1 dimensions with the 
solution 


p{r) 


1 - 

1 + ■ 


(97) 


The energy of this solution is 


E[V]a^0 


271 drr (n'* — p^) 


TT® p 

12 A 




(98) 


The picture is thus just as described in the previous subsection. Since 
the two quark state is accompanied by the vortex of the held V, the 
quarks effectively ‘dig a hole’ in the vacuum. In their immediate vicinity 
the modulus p vanishes, and therefore there is a ‘bag’ of the nonconhning 
state. The radius of this bag is given by the mass of the scalar glueball 
2 ^. 

It is interesting to note that, although for large separation a the energy 
of the string gets contributions from both, the scalar and the pseudoscalar 
glueballs (the modulus and the phase of V), the ‘bag constant’ is deter¬ 
mined solely by the scalar glueball. For small a the phase of V is not 
excited and only p deviates from the vacuum state inside the ‘bag’. This 
is consistent with the common lore that the inside of the bag is distin¬ 
guished from the outside by the value of the condensate. In fact since 
p has vacuum quantum numbers and interpolates in our effective theory 
the scalar glueball, it is naturally associated with the operator which 
has a large overlap with the scalar glueball in QCD. 

The bag we are talking about here arises in a very different situation 
than in the usual bag model j^. There the bag describes the structure 
of the baryon containing light quarks. The radius of the bag in this 
situation is determined by the balance of the vacuum pressure and the 
pressure due to the free motion of the light quarks inside, and in fact 
depends on the quark wave function. In our case the inside of the bag 
contains heavy quarks. Their kinetic energy is small, and we have treated 
them here as static. The radius of the bag thus is determined purely by 
the dynamics of the scalar glueball held and is not sensitive to the state 
of the heavy quarks. This is true for low lying excitations for which the 
radius of the quark state is smaller than the inverse glueball mass. When 
these two scales are comparable presumably the quark pressure will also 
be important and will play a role in the energy balance. Thus in this 
intermediate regime we expect the l/-vortex to be similar to the bag in 
the usual bag model. For states of even larger size the potential between 
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quarks is linear with the fundamental string tension. The bag picture 
should therefore go smoothly into the string picture. 

Three space-time dimensions are unique in many ways. In particular 
they allow massive gauge theories, where the mass of the gauge helds 
does not arise due to the Higgs mechanism. We mean here the Chern- 
Simons gauge theories. The discussion of the magnetic symmetry in 2-1-1 
dimensions would not be complete if we did not discuss this type of models. 
The next section is devoted to such a discussion. 


6 Magnetic symmetry and the Chern- 
Simons term. 

The question of vortices in the Chern-Simons theories is interesting for 
the following reason. The Yang-Mills-Chern-Simons (YMCS) theories are 
completely massive, and so the Wilson loop is expected to have a perimeter 
law behavior |^]. By the argument we have established earlier, it should 
then follow that the magnetic symmetry is unbroken and the spectrum 
should contain vortex states. This is an interesting proposition worth 
investigating. This conclusion would not hold if the magnetic symmetry 
had no local order parameter. If that were the case the symmetry can be 
unbroken, W may have perimeter law and there could still be no magnetic 
vortex states. This is the main question we address in this section. We 
want to examine more carefully the question of locality of the vortex 
operator. 0 

We start our discussion by considering a simpler theory - compact 
electrodynamics with CS term. Conceptual questions here are similar but 
the technical side is much simpler. We find that generically the theory 
does indeed contain local vortex operators and a global discreet magnetic 
symmetry which is unbroken. Nevertheless in the continuum limit there 
are no magnetic vortices in the spectrum. The reason is that the energy 
of such a vortex is logarithmically UV divergent. We find however that 
with a particular scaling of the CS coefficient (logarithmically vanishing 
when UV cutoff is removed) the energy of the vortex becomes Hnite. This 
suggests that the theory may indeed have a phase with finite energy vortex 
states and vanishing photon mass. 

6.1 Compact QED with the Chern Simons term. 

The Lagrangian of Abelian Chern Simons theory in the formal continuum 
limit is 

L = (99) 

The gauge coupling has dimension of mass and k is dimensionless. 
Equations of motion read as dvF^'^ = Kg^e^'^^d^Ap. The mass of the 
gauge particle is M = ng^. The canonical structure of the theory is 
simplest in the Hamiltonian gauge, Aq — 0. 
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The Hamiltonian is 




( 100 ) 


with canonical momenta related to the time derivatives of the fields by 

The gauge helds and canonical momenta form the canonical algebra, 
and the algebra involving the electric fields is 


[Ei{x),Ej{y)] = -iKg eijS {x-y), 
The Gauss law, 


[Ai{x),Ej{y)] = -ig SijS {x-y) 

( 101 ) 


diE^ — ng^ di A j = 0 


( 102 ) 


on a spatial plane generates time independent local gauge transformations. 
The elements of the local gauge group take the form 

17(A) = exp| J (f X X{x)(^diE'— g^ne'^diAj^^ (103) 

such that UAU~^ = A + dX. In the non-compact theory X € R sh ould be 
a single-valued function such that the eigenvalue of the operator (103) on 


physical states is unity. Singular A’s correspond to transformations which 
are in general nontrivial on the physical states. 

Our interest however is in the compact theory. This means that mag¬ 
netic vortices of flux 27r must be physically unobservable. As discussed 
in this amounts to further restricting the physical Hilbert space to 
states which are trivial under the action of the vortex operator. 

In other words certain large gauge transformations must act on the 
physical states trivially in the compact theory as opposed to the non¬ 
compact one. The compact gauge group therefore includes these singular 
gauge transformations in addition to the regular ones, which form the 
gauge group in the noncompact theory. Consider a multi-valued angle 
function 9{x,xo) which is singular at one point and has a discontinuity 
along a straight curve C{xo) that starts at the point xo and goes to in¬ 
finity. The operator of the gauge transformation with this singular gauge 
function creates a magnetic vortex. Its explicit form (after partial inte¬ 
gration and dropping a boundary term owing to the fact that all gauge 
invariant helds decay at inhnity) is 


1 . 


17(a;o)=exp{- J d xdi6{x — xo)\E'' — g Ke^^Aj\} 


(104) 


We have dehned 


di6{x — xo) = di6{x — xo) — 2neijc{x)jS{x — C(xo)) = 


eij{x - xo)j 

(x — xo)^ 


( 105 ) 


c(x)j is a unit vector tangent to the curve C(x, xo)- 
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We may want to include V into the compact gauge group. However 
to be part of the gauge group, it must commute (at least weakly) with 
other elements of the group. One can check explicitly that V(x,xo) does 
not commute with the elements of the noncompact group. To rectify this 
situation we define, following a slightly modified operator 

V{C,xo) = exp^^eij J (f xc{x)jS{x — C{x,xo))Ei{x) (106) 

This operator is merely a “collection” of the electric fields which are per¬ 
pendicular to the curve C{xq). More explicitly one can write it in the 
following form 

Sttz f 

V(C,xo) = exp—eij / dkEj{x) (107) 

9 Jc 

Gauge invariance of this operator, [V, U] = 0, follows immediately. We 
also need to check the commutativity of V{x) with V{y). A straightfor¬ 
ward calculation gives , 


V{Co)V{Ci) 


V{Ci)V{Co)exp 


idtTT K,L(Co, Cl) 


(108) 


Where L{Ci,C 2 ) ~ ±1 if the curves cross each other and L = 0 if they 
don’t. In order that H be a Lorentz scalar the commutator should not 
depend on the curves Ci and C 2 To guarantee this we need to set 

drrK = k £ Z (109) 

We find therefore that the requirement of compactness quantizes the coef¬ 
ficient of the Chern Simons term very much like in the non-Abelian theory 


The commutator 


[B{x), ^"*(* 0)1 = 2Ti7nS^{x - xo)V”^{xo), 


m £ Z. 


( 110 ) 


indicates that (107) creates magnetic vortices of integer strength. Being 
gauge invariant this operator also creates an electric charge 




(fxdiE* = ^ 
47r 




( 111 ) 


Since V has to be included in the gauge group, the magnetic flux and the 
electric charge created by it must be unobservable. Therefore the Hamil- 
ton ian o f the theory must commute with V. The noncompact Hamiltonian 
eq.(lOO) does not quite do the job. It should be modified but in such a way 
that in the continuum limit the same form is recovered for smooth fields. 
The modified Hamiltonian that satisfies these conditions has been sug¬ 
gested in 1^. Since the UV structure is important for our considerations, 
it is most usefully presented in the lattice notations 


Hb = 


1 


a'^g^r 


■E' 


1—Re e 


ina'^ B(x) 


). 


He = 


2 2 

m g 

47r2o2 


E' 


1—Re e 


( 112 ) 


Ij fij E-i {x) \ 

j’ 
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a is the lattice spacing and m,n £ Z and hj is the unit vector parallel to 
the link. The normalization of the electric and magnetic terms is such that 
in the naive continuum limit a —> 0 they reduce to and respectively. 

Using the Gauss’ law one can see that if 2n = fc , the magnetic part 
Hb becomes a combination of the vortex operators V. Therefore without 
loss of generality we assume 2n < k. For m = 1 the electric part of the 
Hamiltonian is also a sum of a fundamental vortex and anti-vortex, and 
we take m > 1. 

Now that we have the formulation of the compact CS QED we can ask 
about the locality properties of vortex operators. The operator V we have 
considered so far is of no interest of itself, since is is trivial on all physical 
states. We thus have to look at the operators which create magnetic flux 
smaller than 27r 

Vp{C,xo) = / dliEj{x), P^Q (113) 

9 Jc 

Here p is a rational number p £ (0,1). The question we are asking is, 
are there such values of p for which 14 > is a gauge invariant local opera¬ 
tor. The gauge invariance with respect to the noncompact gauge group 
is straightforward, since Vp only depends on the electric field, and the 
electric field itself is gauge invariant. However Vp should also commute 
with the ’’fundamental” vortex V, since V is part of the gauge group. 
Therefore we have 


[U, Up] = 0 


kp = I £ Z 


and I < k. 


(114) 


The condition of locality requires that Vp(x) commute with each other at 
different points x and y. This com muta tor also should be independent of 
the contour C in the dehnition eq.(113). 


[Vp{x),Vp{y)] = 0 


kp = r £ Z 


(115) 


Both equations (114) and (115) have to be satisfied for the existence of 
non-trivial local vortices. Whether it is possible or not to satisfy these 
equations clearly depends on the CS coefficient k. For example there are 
no solutions for k — 2 and k = 3 theories. For fc = 4 we can choose I = 2 
and this gives a vortex of vorticity p = 1/2. In general one can solve the 
constraints in the following way. Writing CS coefficient in terms of its 
prime factors, k = gi 92 93 ... Qm, where all qi are not necessarily different, 
one has the following two conditions to satisfy 


qi 92 93--- qmP = l, 


9i 92 93---9mP =r. 


(116) 


The first condition is solved if p divides k, 

p = -i-, where i < m (HI) 

91 92 93--- qi 

Using this in the second condition one can see that the most general form 
of k which allows vortices is 

o 1 

If k = z p — t>2 and t,z£Z (US) 
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For example if fc is a prime number there are no solutions. Generically it 
is easier to find a solution at large values of k. 

The above relations also show that should a solution exist, there is al¬ 
ways a vortex of minimal vorticity. All the other local vortices are simple 
powers of this minimal vortex. For example for k = 36, the above condi¬ 
tions give three solutions (and their integer multiples), p = (1/2, 1/3, 1/6). 
Obviously “p= 1/6” is the minimal vortex. The minimal value of p = ^ 
determines the global magnetic symmetry group of the theory as Z™. 

One last requireme nt th at Vp should satisfy, is locality with respect to 
the energy density eq.(112). In obvious notation 


[hE{x),Vp{y)] = 0, x^y- 


kp 

— G 

m 


(119) 


This can always be satisfied by choosing m = k. To satisfy the other 
condition 

[hB{x),Vp{y)]^O^^GZ (120) 

we can take n = 2. Certainly one can define other Hamiltonians which 
will be compatible with the above conditions. We see therefore that for 
those values of k for which vortex operators are local with respect to each 
other we can always choose the Hamiltonian such that they are also local 
relative to the Hamiltonian density. 

Thus we conclude that for many values of k local physical vortex opera¬ 
tors exist. They are order parameters for a global magnetic symmetry. 
The value of w is determined by k through the solution of the equations 
for minimal p. Thus the argument described above applies and, at least in 
the lattice theory there are vortex states. Calculating their energy in the 
lattice theory is not a simple matter. However the interesting question 
is whether these states survive in the continuum limit. That is to say, 
whether their energy stays finite as the lattice spacing approaches zero. 

In the continuum limit for smooth configurations of the fields the the¬ 
ory is described by the Lagrangian eq. (^^ . However while solving contin¬ 
uum equations we may sometimes encounter field configurations with fast 
variations. For these configurations it is important to take into account 
the compactness of the theory. In particular consider the electric field 
created by the "minimal” vortex operator 


[V{x),Ei{y)] = Ei{y) + ei{x,y), ei{x,y) = —g‘^Kh{y)iS{y-C{x,y)) 

( 121 ) 

where h{y) is the vector tangential the curve C at the point y. Since 
the operator V{x) is local, its only observable action in the compact the¬ 
ory is at the point x. How ever if we just calculate the energy using the 
naive Hamiltonian eq.(lOO) we find infrared divergence proportional to 
the length of the curve C. Clearly if faced with this type of configurations 
in continuum calculations we should subtract this infrared divergence by 
hand. Rather than do this we find it convenient to think about it in the 
following way. Let us split the general electric field configuration into a 
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smooth pie ce an d a piece that contains arbitrary number of strings of the 
type of eq.( 121 ) 

E’' = El^ooth + e* (122) 

and subtract the contribution of a in the Hamilto nian. The only remnant 
of Bi then is in the Gauss’ law, since a of eq.(121) corresponds to a 
pointlike charge ^ at the point x. Thus the smooth field satisfies 

not the naive Gauss’ law, but rather a modified one 


d.El^ootH\^,^^sH^) - = 0 (123) 

In other words we can work entirely in terms of Egmooth if we remember 
that we may allow Gauss’ law to be violated by the presence of pointlike 
charges of charge Kg’^jw. The appearance of w in this way is the only 
remnant of the compactness of the theory. In the following we will work in 
terms of the smooth fields but will drop the subscript smooth for brevity. 

With this caveat in mind, to determine the energy of the magnetic 
vortex in the continuum limit we now should solve the continuum equa¬ 
tions of motion. For a minimum vorticity solution 1/w following | p2| , ^ 
one can take the time independent symmetric ansatz, 

Ai{r) = Bij^[g{r) - —] Ao{r) = h{r) (124) 

w 

The equations of motion read 


g''(r) — —g'(r) — rMh'(r) = 0, h''(r) + -h'(r} — —g'(r) = 0, (125) 

r r r 

where M = g^K. We are looking for the solutions with vorticity 1/w. The 
magnetic field is B = —^g'(r) so we impose (/(O) = 1/w and < 7 ( 00 ) = 0 
and we also demand that the fields decay exponentially at infinity. The 
solution under these conditions can be found as 


g{r) 

A4t 

-- Ki{Mr), 

w 

B{r) 

=- Ko(Mr) 

w 

(126) 

h{r) 

M 

= - Ko(Mr) 

w 

E' = 

- Ki{Mr), 

wr 

(127) 


where Kq,K\ are the Bessel functions.The energy of this vortex follows 
as 


£ = 


TT 


, ,,dh.o 1 ,dn.n^ TTK,,/ , , M\ , 

= (128) 


where A is the ultraviolet cutoff scale. 

This result warrants several comments. First, we see that the energy of 
the vortex is IR finite. This is closely related to the fact that the locality 
of the operator V^/w allowed us to ” vi olate ” the Gauss’ law. Looking at 


the electric and magnetic fields in eqn (127) we indeed see tha t th e naive 
Gauss’ law is violated precisely by the amount allowed by eq.(123). 

Second, the energy is UV divergent. Thus the magnetic vortices do 
not survive in the continuum limit as finite energy excitations. This of 
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course does not contradict our original argument. The Zw symmetry is 
unbroken in the vacuum and excitations carrying the quantum numbers of 
this symmetry are very heavy. The symmetry therefore seems completely 
irrelevant for the low energy dynamics. The curious thing though is that 
although the energy of the vortex is large, it is still much lower than 
the natural ultraviolet scale A. The vortices therefore are not genuine 
ultraviolet objects in the lattice theory, but rather occupy an intermediate 
scale between the UV scale A and the IR scale M. In fact bringing M down 
to zem makes vortices light. At M = ^ln“^ ^ the energy of the vortex is 
hnitetJ. This behavior in fact is very reminiscent of vortices in the Higgs 
phase of the Abelian Higgs model. The mass of such a vortex at weak 
coupling is very large; M oc In where is the mass of 

the massive photon and Mh is the mass of the Higgs particle. However 
as the photon mass decreases, that is as the theory approaches the phase 
transition line, vortices become light. On the phase transition itself they 
in fact become massless and condense in the Coulomb phase. It is not 
unlikely that similar phenomenon occurs in our model. As k decreases 
at fixed the photon becomes lighter, and the mass of the vortex also 
decreases. It could happen that at some value of k the vortices actually 
become massless and drive a phase transition into a phase with broken 
magnetic symmetry. A transition of precisely such type was conjectured 
to happen in the lattice model at fc = 8 in and was seen in the 
variational calculation of pl[|. Of course, within the naive continuum limit 
we consider here we are unable to see such behavior. However the fact 
that the vortices become light within the validity of the naive continuum 
limit is quite suggestive in this respect. Note that if the vortices indeed 
do condense, the magnetic symmetry is spontaneously broken. By 
virtue of the argument given in Sec. 3 this means that the low k phase is 
confining. 


6.2 Non-Abelian theories. 

We now extend this analysis to non-Abelian SU{N) theories, 

L = (129) 

CS coefficient has the well known quantization, 47rK = k 0 . The classical 
equations of motion follow as 


= 0 (130) 


^^The reader may wonder why we are not bothered by the factor k in eq.( |l2q ). After all 
we saw that in the continuum limit naturally > oo. The point is that the flux of the 
minimal vortex generically scales as • 111 ^ oc k . For example ii k = with some integer x, then 
clearly the minimal solution of the eq.( |ll8| ) corresponds to w = x. Thus kluP" is finite in the 
continuum limit. 
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The canonical structure of this theory is similar to the Abelian case. In 
the Ag = 0 gauge 

= + where = Al (131) 

9 ^ 

The Hamiltonian is 

+ where = ie,, (132) 

The canonical algebra is 

[Ei{x),E^j{y)] = -i5‘"*’Kg‘^eij6^{x-y), [A“(f), ^‘’(y)] = -ig^5°-'‘SijS^{x-y) 

(133) 

In terms of the momenta the Gauss law is 

{DiHiY = (134) 

In the non Abelian YMCS theory, the large Wilson loop still commutes 
with the Hamiltonian. This is obvious in the Hamiltonian formalism, since 
the commutation relation between the vector potential and the chromo¬ 
electric field is unaffected by the presence of the CS term. The form of 
the Hamiltonian in terms of Ai and Ei is also the same as without the CS 
term. Since IT is a function of Ai only, its commutator with the Hamil¬ 
tonian is exactly the same as in the theory without CS. Therefore the 
fundamental Wilson loop still generates a symmetry. 

Again our question is whether the theory admits local vortex oper¬ 
ators. In the non-Abelian theory our choices are more limited than in 
compact QED. In the SU(N) Yang Mills theory, the only candidates for 
local operators are those that create quantized flux Q. The vortex 
operator in YM theory is [^, 

Atti f 

V{x)=exp{-^ J dy''eijTr{YEi{y)) (135) 

where the hypercharge generator Y is 

Y = diag(l,l,...,-(Y-1)) (136) 

We have proved in the earlier sections that in SU(N) YM theory, this 
operator despite its nonlocal and gauge non-invariant appearance is in 
fact a local, gauge invariant, Lorentz scalar held Q The way it was 
constructed there was to require that it satishes the ’tHooft algebra 
with the fundamental Wilson loop 

V\x)W{C)V{x) = exp{^n(a;, C')}IT(C) (137) 

with n{x, C) being the linking number on the plane between the point x 
and the closed curve C . 
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One can see that the operator in (135) is also the appropriate vortex 


operator when a Cher n-Sim ons term is included for the gaugefield. The 
commutation relation (137) is still satisfied by the expression (135). How¬ 


ever an additional requirement is that V be gauge invariant, in the sense 
that the matrix elements of V between physical states (those that satisfy 
Gauss’ law) and non-physical states (non-singlet under gauge transforma¬ 
tions) vanish. This means that when we calculate matrix elements of any 
number of operators V We thus have to check that it transforms a physical 
state into another physical state. The difference with the pure YM the¬ 
ory is that the wave function of a gauge invariant state does not depend 
only on the Wilson loops. The physical wave function should satisfy the 
following equation 

i{D.^rnA]^^UjdU^nA] (138) 

The general form of has been determined in 0 in terms of certain 
nonlinear variables. For our purposes we find it more convenient to work 
directly in terms of the vector potentials Ai. Let us take in the form 


= exp{—iS} 

The eikonal S satisfies a linear inhomogeneous equation 

Df^,S[A] = 


(139) 


(140) 


The solution of the homogeneous equation is indeed any functional that 
depends on Wilson loops S'o[W]. We can find a particular solution of 
the inhomogeneous equation using the following argument. 5'[A] must 
be a functional whose change under a standard gauge transformation of 
vector potentials SAi = DiX is proportional to ^ J d^xeijdiA’jX'^. Such a 
functional can be represented as a Chern Simons action on a space with 
a boundary. Let us introduce an additional coordinate r € [—oo, 1] and 
functions of three coordinates Ai{x,T) so that at the boundary r = 1, the 
value of these functions is equal to the value of the vector potentials in 
our theory Ai{x,r = 1) = Ai{x). Let us write the Chern Simons term (in 
the Weyl gauge) on this manifold 


Scs = / dr d^xeijAi{x,T)Aj{x,T) 


(141) 


Under the r independent gauge transformation this action changes by a 
boundary term 


SScs = — J d‘^xeijX'^{x)diAj{x,T = 1) = — J d^xtij\°'{x)diA‘^{x) 

(142) 

which is precis ely of the form required to satisfy eq.(14C). A particular 
solution of eq.(140) is therefore 


K K 

S, = --Scs = -- 


dr / d^xeijAi{x,T)Aj(x,T) 


(143) 
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The introduction of the extra coordinate r and the expression eq.(141) is 
not at all unnatural. One should view this extra coordinate as parame¬ 
terizing a curve in the field space. With this interpretation we have 


dTAj{x,T) = 5Ai 


(144) 


and 


dr 


d^xeijA^{x, T)Aj{x^ r) 


SAi^ij Aj 


(145) 


J —oo J J C 

where the line integral is taken over the trajectory C in the field space 
which ends at the point {Ai{x)}. 

We have thus determined the general form of the wave function of a 
physical state in the YMCS theory to be 


= expli^ 


dr / d xeijA'‘{x,T)Aj{x,T)}'i>o[W] 


«/ — OO «/ 

Now it is straightforward to see how the vortex operator acts on it. 
the action of the vortex operator 


(146) 

Under 


, IL’TT 

V{x) At{y,T)V\x) = Ai{x,T)+—TrYT’^eij 


dzjS‘^{z-C{x,y)) (147) 


Remembering that Ef = —-|- g^KetjA^, we see that the change in 
the phase factor in the wave functional is exactly cancelled by the A- 
dependent term in the vortex operator eq.(^^. 

V{x) V^x) = (148) 


Thus 

V'^[A] = exp{i^ J J (149) 

Clearly a gauge transformed vortex operator Vn has exactly the same 
action on the wave functional 4 '[j4], 

UT = Vh4' (150) 

which establishes gauge invariance of V in the same sense as in the YM 
theory. 

We now can check the locality of the operator V by calculating straight¬ 
forwardly the relevant commutation relation. A simple calculation gives 
\V{x), V{y)] = 0. Thus the operators are local with respect to each other. 
When considering the locality with respect to the Hamiltonian density we 
are faced with the same ambiguity as in the Abelian theory. The electric 
part of naive continuum Hamiltonian is not local relative to V, since Ei is 
shifted by the action of V along the curve C. Just like in the Abelian case 
one should consider a properly regularized version of H in order to be able 
to draw a definite conclusion. In the non-Abelian case such a regularized 
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Hamiltonian is not available. However in the Abelian case we saw that 
there is quite a lot of flexibility in defining such a regularized version. 
In particular we saw that whenever the vortex operators were local with 
respect to each other, we were always able to define the local Hamiltonian 
density. We expect that this situation persists in the non-Abelian theory 
too. 

The situation in the continuum limit is again similar to the Abelian 
case. There are no finite energy solutions of the non-Abelian equations of 
motion which have finite vorticity. The only way to And such solutions 
would be again to relax the Gauss’ law constraint by allowing point like 
charges which correspond to the singular chromoelectric field created by 
V. However again those IR finite configurations will have UV logarithmi¬ 
cally divergent energy. In fact taking Abelian ansatz the YM equations 
of motion reduce to those we considered in the previous section and thus 
lead to the same energy dependence on the UV cutoff. Strictly speaking 
this conclusion is only valid for large enough value of k, since for small k 
quantum corrections to this classical analysis may be large. Thus again it 
is possible that at small k the theory is in a different phase as suggested 
in 0. 


To summarize, we have found that compact CS QED does admit local 
vortex operators for many values of the CS coefficient k. The energy of the 
vortex excitations however generically is logarithmically UV divergent in 
the continuum limit. With a particular scaling of the CS coefficient these 
vortices become light and might condense at small values of k. Our results 
for non-Abelian CSYM theory are similar. Local vortex operators exist, 
but the particles that carry vorticity are heavy in the continuum limit. 

These results are broadly compatible with suggestions made in the 
literature that at low values of the Chern Simons coefficient the YMCS 
theory might undergo a phase transition. If this happens it is very likely 
that this other phase has a broken magnetic symmetry and is therefore 
confining. This is a very interesting possibility which seems worthwhile 
exploring. 


7 Summary 

In this review we have discussed various aspects of the magnetic symmetry 
in 2-1-1 dimensional gauge theories. This symmetry is sometimes called 
topological and is looked at as a somewhat esoteric concept. In fact, as we 
have discussed here, it is not topological at all, in the sense that its charge 
density (or group element per unit volume in the discrete case) is a gauge 
invariant local operator. Moreover in gauge theories without matter fields 
in the fundamental representation it has a local order parameter. It thus 
behaves in all respects as a normal nontopological symmetry. As such 
its mode of realization is rigidly linked to confining properties of the the¬ 
ory. Whenever this symmetry is spontaneously broken, the gauge theory 
is confining, while unbroken realization of magnetic symmetry signifies 
existence of massive gauge bosons. 
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We have constructed the low energy description of gauge theories based 
on the effective Lagrangian for the order parameter of the magnetic sym¬ 
metry. This Lagrangian realises confinement in a very straightforward 
and simple way on the classical level. 

In theories with fundamental matter, the magnetic symmetry is still 
present, but it does not have a local order parameter. As a result it is 
implemented as a local rather than global symmetry. We have shown how 
to implement it in a consistent way in the low energy effective Lagrangian 
in the case of heavy fundamental matter. Interestingly, the resulting ef¬ 
fective Lagrangian provides a bag like description of baryons, as quarks 
surrounded by a ’’bag” of nonconfining vacuum. 

In Chern-Simons theories we proved that the magnetic symmetry is 
also present and is unbroken in the vacuum. Magnetic vortices - the 
excitations carrying the magnetic charge are however infinitely heavy in 
the continuum limit, at least as long as the Chern-Simons coefficient is 
large. For small coefficient there is a possibility that the vortices become 
massless and condense, thus driving the theory into the confining phase. 

In short, magnetic symmetry is a ’’trademark” of 2-1-1 dimensional 
gauge theories, and most qualitative features of their low energy dynam¬ 
ics is determined by the realization of this symmetry. It is hard to over 
emphasise the importance of symmetries, especially in strongly interact¬ 
ing theories, where a symmetry may be our only lead in understanding 
dynamics. Even more so in pure gauge theories, where the magnetic sym¬ 
metry is the only global symmetry of which existence we know. Hopefully 
we will learn how to utilize it to the fullest and will thus gain more infor¬ 
mation about the strongly coupled regime of these theories. 
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